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SECTION  I 
INTRODUCTION 


A.  DISCUSSION 

The  Iranian  Long- Period  Array  (ILPA)  is  located  jar 
Teheran,  Iran.  This  array,  which  became  operational  on  1 May  1976,  is 
a seismic  recording  installation  comprised  o£  a central  recording  station 
and  an  array  of  seven  remote  sites. 

This  report  presents  the  results  of  a continuation  of  the 
evaluation  of  this  installation.  In  the  preliminary  ILPA  evaluation  (Strauss, 
1976),  attention  was  focused  on  data  quality  and  sources  of  data  errors, 
beamforming  gains  in  signal-to- noise  ratio,  first  estimates  of  detection 
capability,  and  M^-rn^  relationships.  Due  to  the  limited  time  available 
between  reception  of  ILPA  data  and  the  end  of  the  contract  period,  no  noise 
analysis  was  performed.  Since  no  recognized  presumed  nuclear  explosions 
occurred  during  May,  1976  (from  which  time  frame  the  data  base  was  drawn), 
it  was  not  possible  to  comment  on  the  ILPA  discrimination  capability. 

In  this  continuation  of  the  ILPA  evaluation,  emphasis  was 
placed  on  improving  and  refining  the  work  started  under  the  preliminary 
evaluation.  With  sufficient  data  now  available,  it  was  also  possible  to 
investigate  the  nature  of  the  noise  field  at  ILPA,  regionalize  the  ILPA 
detection  capability  estimates,  and  consider  the  ILPA  discrimination 
capability. 

B.  THE  SYSTEM 

The  instrumentation  and  operation  of  the  Iranian  Long-Period 
Array  have  been  described  in  detail  in  the  operation  and  maintenance  manual 
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for  the  ILPA  seismic  system  and  in  the  report  on  the  preliminary  evaluation 
of  ILPA  (Strauss,  1976).  Therefore,  it  is  only  necessary  here  to  briefly 
describe  the  system. 

Each  remote  site  of  the  seven-element  array  has  a three-com- 
ponent broadband  seismometer  (KS  36000)  located  in  a 100-meter  deep  bore- 
hole to  reduce  wind- generated  noise.  Each  seismometer  is  a force-balance 
type  which  produces  an  output  proportional  to  earth  accelerations  over  the 
frequency  range  0.  02  - 1.0  Hz.  The  instrumentation  at  each  remote  site  also 
includes  a data  acquisition  subsystem,  a telemetry  subsystem,  and  a power 
subsystem.  The  data  recorded  at  each  site  are  transmitted  directly  to  the 
central  recording  station,  with  one  exception.  Since  site  6 does  not  have 
line-of- sight  with  the  central  recording  station,  the  data  from  site  6 are 
relayed  through  site  5 to  the  central  recording  station.  The  locations  of  the 
remote  sites  are  listed  in  Table  I- 1 and  shown  in  Figure  1-1, 

The  central  recording  station  processes  and  records  data 
received  from  the  seven  remote  sites.  The  instrumentation  housed  in  the 
central  recording  station  includes  the  station  processor,  the  visual  recording 
system,  the  magnetic  digital  tape  recording  system,  and  the  timing,  telemetry, 
and  power  systems.  The  visual  recording  system  converts  the  digital  signals 
from  the  remote  sites  back  to  analog  form  for  display  on  drum  recorders  and 
develocorders.  The  digital  magnetic  tape  recording  system  records  three 
components  of  long-period  motion  from  each  of  the  seven  remote  sites.  This 
system  is  also  used  for  the  tasks  of  providing  data  for  beamforming  and  dis- 
play and  of  editing  data  to  other  tapes. 

The  output  of  the  ILPA  data  recording  system  available  to  this 
evaluation  task  is  the  digital  magnetic  tape  recorded  in  the  satellite  tape  for- 
mat. This  is  a second  digital  magnetic  tape  recording  system  which  records 
three  components  of  long-period  motion  from  each  of  the  remote  sites  and 
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TABLE 
REMOTE  SITE  C 


Site 

Location 

1 

1 (ref) 

35°24'58.  3" 

50°41'19.5" 

2 

35°39'46.  1" 

50°53*51 .5" 

3 

35°28<34.  0" 

51°01'25.5" 

4 

35°14'19.  3" 

50°54'04, 2" 

5 

35°12'46.  2" 

50°34'52.  0" 

6 

35°28'25.  2" 

50°25'32.  2" 

7 

35°42'10.  1" 

50°36'32.0" 

3 


Distance  From  Reference  Site 
(km) 

North  ■ 

East 

0.  0 

0.  0 

27.277 

19.  035 

6.  217 

30.  377 

-19.  536 

19.162 

-22.415 

-9.  830 

5.815 

-23.775 

31.700 

-6.  951 

6 


7 


? 


1 

Reference 

Site 


5 

Relay 


I 

50°4l'l9.  5" 


FIGURE  1-1 

ILPA  SITE  LOCATIONS 
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one  component  of  short-period  vertical  motion  (presently,  from  eito  7)  con- 
tinuously. This  tape  is  800  BPI,  9-track,  recorded  using  two's  complement 
binary  arithmetic.  The  data  are  quantized  at  20.  951  computer  counts  per 
millimicron  (0.  0477  mfi/cc)  of  ground  motion  for  long-period  and  16.  393 
computer  counts  per  millimicron  (0.  061  mjU/cc)  of  ground  motion  for  short- 
period.  The  data  sampling  rate  is  one  sample  per  second  of  long-period  data 
and  twenty  samples  per  second  for  short-period  data. 

C.  EVALUATION  GOALS 

The  specific  goals  of  this  evaluation  are: 

e To  determine  data  quality,  sources  of  data  errors,  and 

reliability  for  each  remote  site  and  for  the  array  as  a whole. 

e To  investigate  the  long-period  noise  field  characteristics  at 

each  remote  site  and  for  the  array  as  a whole. 

• To  investigate  the  signal-to-noise  ratio  gains  due  to  beam- 
forming, 

• To  investigate  site-to-site  signal  similarity. 

• To  determine  estimates  of  the  array  detection  capability  on  a 
regionalized  basis. 

• To  consider  the  discrimination  capability  of  the  array  on  a 
regionalized  basis. 

• To  compare  the  performance  of  ILPA  with  an  SRO  single- site 
station  located  approximately  700  km  east  of  ILPA. 

The  method  of  reaching  these  evaluation  goals  is  as  follows. 
First,  a suitable  suite  of  seismic  events  and  noise  samples  is  created  from 
available  event  lists.  The  procedure  for  selecting  these  events  is  described 
in  Section  II.  This  section  also  describes  the  manner  in  which  the  data  were 
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processed  and  presents  some  estimates  of  data  quality  and  array  reliability 
based  on  data  processing  experience.  The  investigation  of  the  local  noise 
field  characteristics  of  the  array  is  presented  in  Section  III,  The  noise  field 
is  characterized  by  RMS  noise  level,  monthly  RMS  noise  level  trends,  peak 
noise  amplitudes,  and  spectral  content,  The  multiplicity  of  data  provided  by 
seven  sites  also  permits  investigating  the  question  of  propagating  noise  by 
means  of  frequency- wavenumber  spectra.  Section  IV  presents  estimates  of 
signal-to-noise  ratio  gains  due  to  beamforming  and  site-to-site  signal  simi- 
larities. Section  V presents  estimates  of  the  detection  capability  of  the  array 
in  terms  of  the  entire  area  of  interest  and  in  terms  of  specified  sub-regions 
within  this  area  of  interest.  In  the  course  of  estimating  the  array  detection 
capability,  the  effect  of  mixed  events  and  system  malfunctions  on  these  esti- 
mates is  considered.  Section  VI  discusses  the  question  of  discriminating  be- 
tween earthquakes  and  presumed  nuclear  explosions.  The  means  of  discrimi- 
nating used  in  this  section  is  the  Mb-  m^  relationship.  The  final  major  area 
of  this  evaluation  is  found  in  Section  VII,  which  assesses  the  performance  of 
the  array  relative  to  a nearby  single- site  station.  Section  VIII  summarizes 
the  findings  of  this  evaluation  and  suggests  future  avenues  of  investigation. 
Section  IX  lists  the  references  cited  in  this  report.  Finally,  Appendix  A 
describes  the  data  base  used  in  this  evaluation. 
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section  rt 

THE  DATA  BASE 

A.  FORMATION  OF  THE  DATA  BASE 

The  data  base  used  in  thip  continued  evaluation  of  the  Iranian  1 

Long-Period  Array  is  essentially  a subset  of  the  data  base  used  in  the  Seis-  i 

mic  Research  Observatories  evaluation'' (Strauss,  1977).  This  data  base  was 
formed  by  first  selecting  all  events  with  Eurasian  epicenters  as  listed  by  \ 

the  Norwegian  Seismic  Array  (l^ORSAR)  bulletin.  (This  was  the  only  available  ! 

event  list  at  the  time  the  data  'base  was  formed.  ) The  time  frame  selected  was  ■ 

from  22  December  1975  to  30  September  1976,  when  the  NORSAR  bulletin  I 

ceased  to  be  issued.  This  yielded  a list  of  2697  events.  ] 

i 

This  list  was  far  too  large  to  be  successfully  dealt  with  and 

• / [ 
was  therefore  reduced.  First,  the  events  were  grouped  by  seismic  regions 

as  defined  by  Flinn  and  Engdahl  (Flinn  and  Engdahl,  1965).  Those  regions 

containing  fewer  than  80  events  were  dropped  from  the  event  list.  Next,  the 

remaining  seismic  regions  were  broken  into  their  geographic  sub- regions, 

again  as  defined  by  Flinn  and  Engdahl.  Those  sub-regions  which  contained 

only  a few  events  and  which  were  separated  from  the  main  body  of  event-  j 

containing  sub- regions  in  the  region  were  then  rejected  and  their  contained  | 

events  deleted  from  the  event  list.  The  remaining  regions  (shown  in  Figure  j 

II- 1)  were  used  in  this  evaluation.  The  regions  have  been  renumbered  but  j 

have  the  same  boundaries  as  those  of  Flinn  and  Engdahl,  with  the  exception  j 1 

of  their  region  19,  which  was  sub-divided  to  form  regions  1 and  2 of  Figure  j 

II- 1.  The  subregions  composing  each  region  are  listed  in  Table  II-l. 
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At  this  point,  the  data  base  contained  approximately  2300 
events.  To  bring  this  down  to  a manageable  size,  events  were  systematically 
deleted  (every  other  event,  two  out  of  every  three  events,  as  needed)  from 
those  regions  containing  more  than  120  events  until  they  each  contained  approx- 
imately 120  events.  This  reduced  the  data  base  to  1472  events. 

The  1LPA  data  base  was  created  from  this  data  base  by  selec- 
ting all  events  for  which  the  array  was  nominally  operational.  This  resulted 
in  a data  base  of  497  events  occurring  between  1 May  1976  and  30  September 
1976.  Finally,  this  data  base  was  increased  to  613  events  by  selecting  events 
from  the  above-defined  regions  as  listed  by  the  National  Earthquake  Infor- 
mation Service  (NEIS)  event  bulletin  from  1 October  to  20  December  1976. 

The  event  parameters  of  this  data  base  are  presented  in  Appendix  A. 

In  the  preliminary  1LPA  evaluation,  the  data  base  consisted 
of  281  events  from  May  1976.  For  the  current  data  base,  these  events  were 
subjected  to  the  same  regionalizing  procedure  as  the  new  events,  reducing 
the  total  for  May  1976  from  281  to  213.  The  67  events  of  the  old  data  base 
which  were  not  included  in  the  new  data  base  were  rejected  because  their 
epicenters  were  not  in  regions  used  in  this  evaluation.  These  67  events 
will  be  used  in  the  discussion  of  data  quality  and  array  reliability  however, 
since  these  matters  are  independent  of  region. 

Formation  of  the  noise  data  base  will  be  discussed  in  detail 
in  the  section  on  noise  analysis  (Section  III),  In  brief,  the  noise  data  base 
was  formed  by  searching  the  NORSAR  and  NEIS  event  bulletins  for  daily 
time  intervals  of  at  least  one  hour  in  duration  during  which  no  seismic  sig- 
nals would  be  expected  to  arrive  at  the  array.  Noise  edits  were  then  cre- 
ated by  selecting  data  segments  4096  seconds  in  length  sampled  at  a rate  of 
one  sample  per  two  seconds  from  these  intervals. 
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B.  DATA  PROCESSING  METHOD 

The  computer  processing  of  seismic  signals  and  noise  samples 
can  be  divided  into  two  functional  stages  - the  basic  processing  stage  and  the 
analysis  processing  stage. 

The  basic  processing  stage  is  outlined  in  Figure  II-2.  The  de- 
sired time  segments  of  signal  or  noise  are  first  edited  from  the  merge  tapes 
created  by  Teledyne/Geotech.  (These  merge  tapes  each  contain  three  consec- 
utive days  of  1LPA  data  recorded  at  1600  BPI.)  The  output  from  the  edit  pro- 
gram consists  of  trace  mean  information  for  each  component  of  each  site, 
quality  control  information,  and  a permanent  hold  tape  of  the  edited  data.  The 
quality  control  information  consists  of  messages  indicating  the  presence  of 
parity  errors,  timing  errors,  clipped  data,  and  spiked  data.  Also  printed  out 
is  a summary  of  segment  powers,  which  can  be  used  to  determine  bad  sites. 
(Bad  sites  are  those  sites  which  are  dead,  contain  uncorrectable  spikes,  or 
display  abnormally  high  or  low  power  levels.  ) At  this  point,  the  analyst  uses 
the  quality  control  information  to  guide  him  in  selecting  the  array  sites  to  be 
input  to  the  beamforming  program. 

The  next  major  program  of  the  basic  processing  stage  performs 
trace  mean  removal,  rotation  of  the  data  from  thoir  recorded  vertical,  north, 
east  (V,  N,  E)  configuration  to  a vertical,  transverse,  radial  (V,  T,  R)  configu- 
ration, and  beamforming  of  the  good  sites.  Rotation  of  the  data  separates  the 
surface  waves  recorded  on  the  horizontal  components,  theoretically  resulting 
in  two  components  of  Rayleigh-wave  motion  (V  and  R)  and  one  component  of 
Love-wave  motion  (T).  Noise  samples  retain  their  V,N,E  configuration.  Both 
the  edit  and  beamforming  programs  operate  on  one  128-point  data  segment  at 
a time,  continuing  until  the  desired  data  length  has  been  processed. 

The  beamforming  program  operates  by  computing  arrival 
time  delays  at  each  site  relative  to  a reference  site  (site  1)  using  fixed 
velocities  of  4.  0 km/ sec  for  Love  waves  and  3.  5 km/ sec  for  Rayleigh  waves. 
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The  data  from  each  site  are  then  time- aligned  using  these  time  delays, 
summed,  and  scaled.  This  process  produces  a signal-to-noise  ratio  im- 
provement by  suppressing  random  noise  by  a factor  approximately  equal  to 
the  square  root  of  the  number  of  sites  used,  ThiB  program  outputs  three 
components  of  motion  from  a reference  site  (usually  site  1 ) and  three  beams 
(V,  T,  R)  to  a permanent  hold  tape  with  appropriate  annotation. 

The  final  program  of  the  basic  processing  stage  performs 
bandpass  filtering  (0.  023  - 0.  059  Hz  passband)  of  the  reference  site  and 
beam  traces  and  outputs  plots  suitable  for  analysis. 

The  various  programs  used  in  the -analysis  processing  stage 
will  be  described  in  the  sections  on  noise  analysis  and  signal  analysis. 

C.  ARRAY  RELIABILITY 

The  first  point  to  be  considered  in  any  assessment  of  array 
reliability  is  the  percentage  of  time  for  which  the  array  is  nominally  oper- 
ational. Table  II-2  shows  that  this  factor  varies  greatly  from  month  to  month. 
(These  percentages  represent  the  percentage  of  each  month  for  which  data 
was  received  at  the  Seismic  Data  Analysis  Center  in  Alexandria,  Virginia. ) 

In  the  first  thirteen  months  of  operation,  ILPA  suffered  from  three  major 
problems.  Between  22  June  1976  and  6 August  1976,  the  array  was  down  due 
to  a malfunction  of  the  air  conditioner  compressor,  which  allowed  temper- 
atures to  rise  above  the  operational  limits  of  the  CRS  hardware.  The 
second  down  period,  from  18  September  to  1 November  1976,  was  due  to 
problems  with  the  satellite-format  tape  recorders.  (From  the  authors' 
point  of  view,  the  array  is  down  when  no  data  is  available  at  the  Seismic 
Data  Analysis  Center  in  Alexandria,  Virginia,  regardless  of  whether  the 
remainder  of  the  sensor  - transmitter  - recorder  system  comprising  the 
array  is  operational.  Thus,  the  array  is  considered  to  be  down  when  the 
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TABLE  II- 2 

PERCENTAGE  OF  TIME  ILPA  IS  OPERATIONAL 


May 

1976 

100% 

June 

1976 

*«  7% 

July 

1976 

3% 

Aug. 

1976 

87% 

Sept. 

1976 

50% 

Oct. 

1976 

0% 

Nov. 

1976 

100% 

Dec. 

1976 

68% 

Jan. 

1977 

0% 

Feb. 

1977 

79% 

Mar. 

1977 

87% 

Apr. 

1977 

100% 

May 

1977 

81% 

Average 

vO 

11 

| . satellite-format  tape  recorders  are  down. ) The  third  down  period,  from 

! 20  December  1976  to  2 February  1977,  was  due  to  problems  with  the  thermal 

| electric  generators  and  the  satellite -format  tape  recorders. 

Overall,  in  the  thirteen  month  period  shown  in  Table  II- 2, 

| , the  average  up  time  percentage  was  only  647c.  However,  as  hardware  prob- 

lems of  the  types  described  above  are  found  and  corrected,  one  can  expect 
| j this  value  to  increase.  An  indication  that  this  is  the  case  may  be  found  in 

Table  II-2,  where  the  average  up  time  for  the  last  four  months  is  87%. 

1 i 

\ ' Even  when  the  array  functions  as  designed  and  tapes  in  the 

i I i satellite  format  are  received  at  the  Seismic  Data  Analysis  Center,  two  fac- 

| 1 - tors  may  prevent  the  analyst  from  processing  a desired  time  frame.  First, 

| j , no  data  is  recorded  during  the  interval  between  the  time  one  tape  is  filled 

| ' ' and  the  next  is  mounted  on  the  tape  drive.  (The  average  time  gap  between 

| j | tapes  is  ten  minutes.  ) Second,  a tape  is  occasionally  received  which  is  un- 

f ' 1 readable  due  to  such  factors  as  poor  tape  quality  or  dust  and  dirt  on  the  tape, 

j When  either  of  these  occurs,  the  analyst  cannot  examine  the  time  frame 

1 in  which  he  is  interested. 


i t 


n 


Out  of  680  event b examined  (613  of  the  current  data  base 
plus  67  additional  events  from  the  data  base  of  last  year),  87  could  not  be 
processed  due  to  gaps  in  the  recorded  data  or  unreadable  data.  Thus,  it 
appears  that  even  if  the  array  is  nominally  operational  100  percent  of  the 
time,  the  analyst  can  expect  to  be  unable  to  process  approximately  12 
percent  of  the  data.  However,  it  should  be  noted  that  45  of  the  87  events 
not  processed  were  lost  to  unreadable  tapes,  a problem  which  should  not  be 
too  difficult  to  rectify.  Correction  of  this  problem  would  bring  the  loss  rate 
down  to  approximately  6 percent. 

Even  after  an  event  has  been  successfully  processed,  system 
malfunctions  can  prevent  the  analyst  from  determining  whether  the  data  shows 
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the  desired  signal,  seismic  noise,  or  a mixed  event.  These  system  mal- 
functions take  the  form  of  spikes,  glitches,  and  data  drop-outs.  These  sys- 
tem malfunctions  degraded  28  of  the  680  events  examined  last  year  and  this 
year  for  a loss  rate  of  4 percent. 

The  overall  ILPA  reliability  estimate  is  determined  from  the 
above  causes  of  loss  of  data.  This  estimate  is  expressed  as: 

Reliability  = Probability  (array  is  operational)  * Probability 
(data  is  recorded  and  readable)  # Probability 
(no  malfunctions). 

The  worst-case  reliability  assumes  no  improvements  in  the  up-time  or 
quality  of  data  recording.  This  estimate  is: 

Reliability  = (0.  64)  * (0.  88)  * (0.  96)  = 0.  54. 

If  one  assumes  the  improvement  in  operational  reliability  sug- 
gested by  the  last  four  months  shown  in  Table  II- 2 is  permanent,  the  array 
reliability  becomes: 

Reliability  = (0.  87)  * (0.  88)  * (0.  96)  = 0.  73. 

Finally,  there  is  no  reason  why  all  recorded  data  should  not 
be  readable.  If  proper  tape  handling  procedures  in  recording  and  shipment 
are  observed,  the  array  reliability  can  be  estimated  as: 

Reliability  = (0.  87)  * (0.  94)  * (0.  96)  = 0.  79. 

This  is  the  best- case  estimate. 

D.  SITE  REJECTION  STATISTICS 

Another  indication  of  the  performance  quality  of  ILPA  is  the 
number  of  sites  considered  to  be  acceptable  for  inclusion  in  beamforming. 

A site  is  rejected  from  the  beamforming  process  if  any  of  the  following  oc- 
curs on  one  or  more  site  components: 


• The  component  is  dead 

• The  data  contain  uncorrectable  spikes  or  clipped  data 

• The  data  contain  power  surges  raising  the  128-point  processing 
segment  powers  more  than  ten  times  the  power  of  the  preceding 
and  following  segments 

• The  data  contain  128-point  processing  segments  with  aero  pow- 
er (data  drop-outs) 

s The  data  contain  128-point  processing  segment  powers  consist- 

ently higher  than  the  segment  powers  of  corresponding  com- 
ponents at  other  sites. 

The  decision  as  to  whether  to  accept  or  reject  a site  is  made  using  the  edit 
printout,  which  lists  segment  powers  for  each  component  of  each  site  of  the 
array. 

Of  the  680  events  in  the  combined  data  bases  of  last  year  and 
this  year,  unreadable  or  missing  data  prevented  the  processing  of  87  events. 
The  number  of  times  each  site  was  rejected  from  beamforming  for  the  593 
processed  events  is  summarized  in  Table  11-3.  (Note  that  site  1 did  not  be- 
come operational  until  5 May  1976  and  site  3 did  not  become  operational  until 
8 September  1976.  ) The  data  of  Table  II- 3 show  that  site  6 had  the  highest 
rejection  rate.  In  almost  all  cases,  the  reason  for  rejecting  this  site  was 
transmission  emrs  occurring  during  the  relaying  of  data  from  site  6 through 
site  5 to  the  central  recording  station.  These  transmission  errors  were  fair- 
ly uniformly  distributed  throughout  the  time  period  covered  by  the  data  base. 
The  transmission  errors  manifested  themselves  in  the  data  as  power  surges, 
spikes,  and  data  drop-outs. 

Table  U-4  presents  the  reasons  for  rejection  of  each  site. 

Single  segment  power  surges  and  entire  component  power  surges  are  pre- 
sented jointly  as  'POWER  SURGES'.  The  heading  'DEAD  COMPONENT' 


TABLE  II- 4 

REASON  FOR  SITE  DELETION 


Percentage  of  Total  Number  of  Rejectlona 


Reason 

1 

2 

3 

4 

5 

6 

7 

Power  Surge 

76% 

36% 

50% 

94% 

34% 

49% 

36% 

Dead  Component 

11% 

57% 

50% 

5% 

61% 

39% 

49% 

Spikes  /Clipped  Data 

13% 

7% 

0% 

1% 

5% 

12% 

15% 

3 


includes  data  drop-outs  and  dead  components.  Following  site  6,  site  ^'/hiad 

• . • v#  fl 

the  highest  rejection  rate.  Table  H-4  indicates  that  the  primary  reason  for 
rejecting  this  site  was  power  surges.  Sites  3 and  5 had  the  next  highest  re- 
jection rates,  For  these  sites,  power  surges  and  dead  component*  caused 
about  the  same  number  of  rejections.  Sites  1,  2,  and  7 had  th&  lowest  re- 
jection rates.  Power  surges, accounted  for  the  majority  of  rejections  of  site 
1.  At  sites  2 and  7,  dead  components  accounted  for  slightly  more  rejections 
than  did  power  surges.  In  all  cases,  spikes  and  clipped  data  accounted  for 
only  a small  number  of  site  rejections. 

E.  PROBABILITY  OF  MIXED  EVENTS  ./ 

He  manner  in  which  the  data  base  for  the  evaluation  of  the 
Iranian  Long  Period  Array  was  selected  produces  a data  base  which  is  essen- 
tially a random  sampling  of  Eurasian,  seismic  events.  Since  no  effort  was 
made  to  exclude  obviously  mixed  events,  the  mixed  event  statistics  derived 
from  analysis  of  this  data  base  should  accurately  reflect  the  frequency  of 

occurrence  of  mixed  event**  (,/k'mixed  event  is  any  event  whose  waveform* 

'v 

are  masked  or  interfered  with  by  the  waveforms  of  a second  event.  ) 

Table  II- 5 p^sents  the  mixed  event  statistics  and  derived 
probability  of  mixing  fo.r  f'ovir  bodywave  magnitude  (m^)  ranges.  These  values 
are  compiled  for  ranges  in  rather  than  for  individual  values  to  provide 
more  robust  statistics.  (It  seems  reasonable  to  assume  that  there  will  be  no 
major  changes  hr*  the  probability  of  mixing  from  one  m^  value  to  the  next.  ). 
This  table  indicates  that,  as  one  might  expect,  the  probability  of  mixing  very 
gradually  dr'ops  as  the  bodywave  magnitude  increases.  This  gradual  decrease 
in  the  probability  of  mixing  is  due  to  beamforming  the  data,  which  suppresses 
off-a*iijn';uth  signals  so  that  events  which  are  mixed  on  the  reference  site  can 
sometimes  be  detected  on  the  beam  data. 
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TABLE  II- 5 

ILPA  BEAM  DATA  MIXED  EVENT  STATISTICS 


Range 

■■  "■  • — —— — 

Number  of 
Mixed  Event* 

Number  of 
Non-Mixed  Event* 

Probability  of 
Mixed  Event 

3.  1-3.  5 

27 

63 

0.  30 

3.6-4.  0 

52 

114 

0.  31 

4. 1-4. 5 

20 

85 

0.  19 

4.6-5.  0 

12 

66 
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SECTION  III 
NOISE  ANALYSIS 


A.  DISCUSSION 

The  goal  of  this  section  is  to  determine  the  long-period  noise 
field  characteristics  of  the  Iranian  Long  Period  Array.  Both  single- site  and 
beamformed  data  will  be  examined  in  order  to  estimate  the  effects  of  the  beam- 
forming process  on  the  noise  field.  The  major  effort  will  be  concentrated  on 
the  noise  in  the  0.  023-0,  059  Hz  signal  window.  This  window  is  used  in  the 
noise  analysis  to  permit  one  to  consider  the  effect  of  noise  level  on  detection 
capability.  Some  information  on  the  noise  field  at  frequencies  outside  this 
window  is  presented  by  average  noise  amplitude  spectra. 

The  analysis  of  the  noise  field  is  divided  into  two  parts  termed 
the  basic  noise  analysis  and  the  extended  noise  analysis.  The  basic  noise 
analysis  covers  the  points  common  to  the  analysis  of  the  noise  field  at  a single- 
site station  such  as  one  of  the  Seismic  Research  Observatories.  These  points 
are  RMS  noise  amplitudes,  RMS  noise  trends,  spectral  content  of  the  noise, 
and  peak  noise  amplitudes.  The  extended  noise  analysis  covers  those  points 
which  can  only  be  investigated  with  array  data.  These  points  are  multichannel 
noise  coherence  and  propagating  noise. 

The  noise  sample  data  base  was  formed  by  searching  the 
NORSAR  event  lists  for  daily  time  intervals  of  at  least  one  hour  in  duration 
during  which  no  seismic  signals  would  be  expected  to  arrive  at  the  array. 

Data  segments  4096  seconds  in  length  were  then  selected  from  these  intervals 
and  processed  as  described  in  Section  H.  After  processing,  the  data  were 
plotted  and  visually  examined  for  unreported  signals  and  system  malfunctions. 
If  either  of  these  was  found,  the  noise  sample  was  rejected  unless  the  signal 
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or  malice'  n occurred  near  tho  start  or  end  of  the  sample.  In  this  case, 

2048  seconds  of  the  noise  sample  were  ufied. 

The  definition  of  seismic  noise  used  in  this  report  is  most 

clearly  stated  by  Enders  A.  Robinson  (Robinson,  1 967)* 

"Any  ground  motion  that  is  not  caused  by  an  explosion  or  an 
earthquake  is  usually  regarded  as  ambient  seismic  noise. 

The  predominant  components  of  such  seismic  noise  are  sur- 
face-generated microseisms  that  originate  from  meteorlogical, 
hydrodynamic,  or  cultural  sources,  Such  microseisms  chiefly 
propagate  along  the  surface  of  the  earth  as  Rayleigh  waves.  " 

B.  BASIC  NOISE  ANALYSIS 

The  goals  of  the  basic  noise  analysis  were  to  estimate  the  long- 
period  RMS  noise  levels,  the  peak  25- second  noise  amplitudes,  and  the  spec- 
tral content  of  the  noise  for  each  of  the  three  components  of  motion  (vertical, 
north,  and  east)  as  recorded  at  each  array  site  and  for  the  beamformed  noise 
sample.  The  overall  evaluation  time  frame  was  from  1 May  1976  to  29  April 
1977.  This  was  used  to  study  the  reference  site  and  beam  noise.  Thirty  days 
of  noise  data  were  selected  from  within  this  time  frame  to  study  the  noise  at 
the  individual  sites. 

After  the  preliminary  processing  and  visual  inspection,  the 
noise  samples  were  input  to  a program  which  performs  the  following  functions} 

i 

e Compute  RMS  noise  values  uncorrected  for  instrument  re- 

sponse in  the  0.  023-0.  059  Hz  passband. 

e Measure  zero-to-peak  25- second  noise  amplitudes. 

e Compute  the  power  spectrum  for  each  component  of  motion  of 

the  noise  sample  and  smooth  to  128  frequencies. 

e Compute  RMS  noise  values  corrected  for  instrument  response 
in  the  0.  023-0.  059  Hz  passband, 
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The  program  output  all  measured  and  computed  values  on 
J punched  cards  to  facilitate  data  input  to  succeeding  plot  programs.  The 

25- second  noise  amplitudes  were  measured  as  the  maximum  25-second 
noise  amplitude  of  each  component  of  motion  of  the  noise  sample. 

The  mean  long-period  RMS  noise  values  in  millimicrons  as 
measured  in  the  0.  023  - 0.  059  Hz  passband  are  presented  in  Table  111-1. 

These  values  are  uncorrected  for  instrument  response.  The  individual 
RMS  values  were  not  plotted  due  to  the  relatively  small  number  of  samples 
used.  In  every  case  except  one  (site  1 north)  the  vertical  component  displays 
somewhat  lower  mean  RMS  noise  values  than  do  the  horizontal  components. 

| There  appears  to  be  no  correlation  between  the  mean  noise  levels  at  the 

sites  and  the  potential  cultural  noise  sources  described  in  the  final  report 
j on  the  installation  of  ILPA  (Texas  Instruments,  1977).  For  example,  an 

active  manganese  mine  is  described  as  being  located  3 km  from  site  3, 

| which  has  mean  RMS  noise  values  at  about  the  median  for  the  seven  sites. 

It  would  appear  that  if  the  types  of  cultural  activity  in  the  area  of  the  array 

j produce  long-period  transient  noise  trains,  these  noise  trains  do  not 

materially  affect  the  mean  RMS  noise  levels. 

1 The  long-period  RMS  noise  values  in  millimicrons  for  the 

, . vertical,  north,  and  east  components  of  motion  are  plotted  versus  Julian 

I I day  for  the  reference  site  (site  1)  in  Figure  III- 1 and  for  the  beam  in  Figure 

I QI-2.  In  these  figures,  RMS  noise  values  for  consecutive  days  are  connected 

) by  solid  lines.  Gaps  of  one  or  two  days  in  length  are  due  to  bad  noise  samples. 

, ■ The  three  large  gaps  are  due  to  array  system  failures,  the  causes  of  which 

l were  discussed  in  Section  II.  The  monthly  RMS  noise  level  trends  which 

j were  derived  from  these  RMS  noise  values  are  shown  in  Figure  lH-3for 

the  reference  site  data  and  Figure  IU-4  for  the  beam  data.  Due  to  array 
( down  time,  it  is  difficult  to  determine  how  the  noise  level  varies  with  time. 

It  appears  that  the  noise  levels  rise  in  the  winter  months  and  decrease  with 

( 
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FIGURE  m-1 

REFERENCE  SITE  (SITE  1)  RMS  NOISE  IN  MILLIMICRONS 


BEAM  RMS  NOISE  IN  MILLIMICRONS 


RI1S  AtfLITUDE  IN  HU 


I I 

( j FIGURE  III- 3 

1 1 REFERENCE  RITE  (SITE  I)  RMS  NOISE  LEVEL  TRENDS 

I! 


m-7 


the  onset  of  spring.  The  only  significant  difference  between  the  reference 
site  data  and  beam  data  as  portrayed  by  Figures  UL-1  to  III-4  is  the  absolute 
noise  level.  Table  XII-2  compares  the  overall  reference  Bite  and  beam 
mean  RMS  noise  levels.  This  table  shows  that  beamforming  reduces  the  RMS 
noise  level  in  the  0.  023  - 0.  059  Hz  passband  by  approximately  6. 1 dB  bn  the 
vertical  component,  4.  8 dB  on  the  north  component,  and  3.  6 dB  on  the  east 
component. 

Table  III- 3 presents  the  statistics  for  the  25-second  noise 
amplitudes  measured  on  each  noise  sample.  The  values  used  tb  compute 
these  statistics  were  the  largest  25-second  noise  amplitudes  of  each  noise 
sample  measured  from  zero-to-peak  in  millimicrons.  The  means  and 
standard  deviations  in  this  table  are  for  the  base  ten  logarithms  of  these 
amplitudes.  These  values  present  another  measure  of  noise  suppression 
due  to  beamforming.  The  peak  25- second  noise  amplitudes  are  reduced  by 
approximately  5.  6 dB  on  the  vertical  component,  5.  2 dB  on  the  north  com- 
ponent, and  5.  0 dB  on  the  east  component  due  to  beamforming.  These 
values  compare  fairly  well  with  the  equivalent  values  for  RMS  noise  suppression 
on  the  vertical  and  north  components.  The  lower  RMS  noise  suppression  for  the 
east  component  (3.  6 dB)  suggests  that  some  range  of  frequencies  on  the  east 
component  in  the  0.  023  - 0,  05  9 Hz  passband  contains  noise  energy  which  is  to 
some  extent  correlated  from  site  to  site  and  which  therefore  is  not  as  well 
suppressed  by  beamforming  as  noise  at  other  frequencies. 

Figure  III- 5 shows  average  RMS  noise  amplitude  spectra  for 
each  of  the  remote  sites  of  the  array  and  for  the  beamformed  noise  data. 

These  spectra  are  uncorrected  for  instrument  response.  They  were  computed 
by  averaging  the  spectra  measured  from  each  noise  sample  and  converting 
these  average  spectra  to  RMS  amplitudes  using  Parsevai's  formula 
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TABLE  III- 2 

COMPARISON  OF  SINGLE-SITE  AND  BEAM  RMS  ttOISE  AMPLITUDES  j I 

(UNCORRECTED  FOR  INSTRUMENT  RESPONSE)  l-J 

(ALL  VALUES  ARE  IN  MILLIMICRONS) 

I! 


V 

N 

E 

L 

No.  of 

Type 

Mean 

S.D.  # 

Mean 

S.  D.  # 

Mean 

S.  D.* 

Samples 

Ref.  Site 

9.  58 

2.  19 

9.  34 

2.  56 

9.  86 

2.98 

9T 

Beam 

4.  77 

1.44 

5.  36 

1.  82 

6.  51 

2.31 

113 

* Standard  Deviation 


HI- 10 


0 

D 


.TABLE  HI- 3 

MAXIMUM  2 5 -SECOND  NOISE  AMPLITUDES 

[ZERO-PEAK*] 
AMPLITUDE  J 


Mean 

Standard 

Deviation 

V 

1.44 

0.13 

Reference  Site 

N 

1.45 

0.14 

* 

E 

1.48 

0.  15 

V 

1. 16 

0.  15 

Beam 

N 

1.  19 

0.  14 

E 

1.  23 

0.  15 

Suppression  of  25-Second  Noise  Amplitudes  by  Beamforming: 
V:  5.  60  dB  N:  5.  20  dB  E:  5.  00  dB 

* Zero-to-Peak  amplitudes  were  measured  in  millimicrons. 
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where  Ai  - 

I A(fi)l2  = 

C^) 


b a 


the  elemental  frequency  interval  0f*O.  001953  He), 

the  discrete  Fourier  transform  spectral  density  es- 
timate at  frequency  f^, 

the  instrument  response  correction  at  frequency  f^, 
the  initial  frequency  index,  and 
the  final  frequency  index. 


Since  the  RMS  amplitude  at  each  discrete  frequency  was  desired,  in  this 
case  a * b.  Also,  since  no  instrument  response  corrections  were  made, 
C(f^)  * 1 for  all  frequencies. 


For  all  sites  except  site  1,  the  RMS  noise  amplitudes  for  per- 
iods greater  than  25  seconds  are  higher  on  the  horizontal  components  than  on 
the  vertical  component.  If  the  noise  on  the  horizontal  components  at  these 
periods  contains  a relatively  large  amount  of  coherent  energy  in  comparison 
to  the  vertical,  the  difference  between  the  reduction  in  noise  due  to  beam- 
forming for  the  vertical  and  horizontal  components  will  be  explained. 


Figures  III- 6 and  HI-7  show  the  average  RMS  amplitude 
spectra  of  all  noise  samples  for  the  reference  site  and  beam  data  respect- 
ively. No  instrument  response  corrections  were  applied  to  these  data.  The 
left-hand  side  of  each  figure  shows  the  average  RMS  amplitude  spectra  for 
the  three  components  of  motion.  The  right-hand  side  of  each  figure  showe  the 
log  RMS  amplitude  spectra.  The  vertical  bars  on  these  spectra  represent 
plua-or-minus  one  standard  deviation  of  the  noise.  These  bars  represent  the 
day-to-day  variability  of  the  noise  at  each  period.  The  following  points  should 
be  noted  from  these  figures: 


■J 

■I 


e The  smallest  decreases  in  the  amplitude  spectra  due  to  beam- 
forming occur  at  the  shorter  periods.  The  most  noticeable 
effect  of  beamforming  at  periods  less  than  25  seconds  is  to 
make  the  vertical  and  horizontal  spectra  nearly  identical. 
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• At  periods  greater  than  25  seconds,  beamforming  lowers  the 
vertical- component  spectrum  much  more  than  the  horizontal 
spectra.  The  east- component  spectrum  shows  relatively  little 
change  due  to  beamforming. 

• The  reference  site  shows  greatest  day-to-day  variation  in  the 
noise  at  periods  greater  than  approximately  28  seconds. 

• Beamforming  greatly  decreases  the  day-to-day  variation  in  the 
noiBe  at  periods  between  14  and  28  seconds, 

• Beamforming  has  very  little  effect  on  the  day-to-day  variation 
in  the  noise  at  periods  greater  than  28  seconds. 

A general  conclusion  which  may  be  drawn  from  the  preceding 
observations  is  that  the  signal-to-noise  ratio  may  be  greatly  enhanced  by 
changing  the  bandpass  filter  limits  from  0.  023  - 0.  059  Hz  (16.  9 - 43.  5 sec- 
onds period)  to  0.  033  - 0.  050  Hz  (20  - 30  seconds  period).  This  will  sharply 
curtail  the  effects  of  the  noise  at  the  microseismic  peak  and  at  periods  beyond 
30  seconds  on  the  signal-to-noise  ratio.  Since  the  majority  of  the  signal 
energy  lies  within  the  0.  033  - 0.  050  Hz  pas  stand,  the  signal  should  not 
be  materially  affected  by  this  change  in  fiLter  limits. 


C.  EXTENDED  NOISE  ANALYSIS 

This  subsection  deals  with  the  types  of  noise  analysis  which 
can  only  be  performed  with  array  data,  i.  e.  , coherence  of  the  noise  and 
directionality  of  the  noise. 

Two  main  programs  were  used  to  carry  out  this  phase  of 
the  noise  analysis.  The  first  program  computes  crosspower  spectral 
matrices  from  the  time-domain  noise  data.  This  is  performed  in  the  follow- 
ing manner.  The  edited  data  In  its  recorded  vertical,  north,  east  configuration 
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is  entered  into  the  program  one  128-point  (256  second) segment  at  a time. 

The  program  then  removes  the  trace  mean  calculated  in  the  edit  program 
from  the  data  segment  and  Fourier  transforms  the  data  segment.  To  pro- 
vide greater  frequency-to- frequency  stability,  a three-point  Hanning  function 
is  applied  to  the  transformed  data,  ThiB  Hanning  function  may  be  expressed 
as 

t;  (f)  = +\Timl(£)  + iTi(f)  + *Ti+1(f) 

where  T^(f)  is  the  input  transformed  data  at  frequency  f 
T^(f)  is  the  output  transformed  data  at  frequency  f 
i = frequency  index, 

At  this  point,  the  crosspower  spectral  matrices  are  computed. 
One  matrix  is  computed  for  each  frequency  within  the  specified  passband 
for  each  component  and  site,  The  program  calculates  these  crosspower 
spectral  matrices  [$]  for  the  desired  sites  i,  j from  the  complex  transformed 

■A 

data  X at  each  frequency  f.  Looping  on  the  number  of  transformed  segments 
(NSEG)  stacks  the  matrix  over  the  entire  data  trace.  In  equation  form 

NSEG 

V>*£  xi«V  • 

Each  element  of  the  matrix  is  scaled  to  account  for  the  number 
of  transform  segments  over  which  the  matrix  was  accumulated  and  to  convert 
from  computer  counts  to  millimicrons.  The  scale  factor  is: 

SCALE  3 Tnptshhsegmq) 

where:  fit  = the  sampling  interval 

NPTS  - thii  nurnoer  of  points  in  each  segment 
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NSEG  si  the  number  of  segments  used,  and 

Q s:  the  quantization  in  computer  counts  per  millimicron. 

Finally,  the  scaled  crosepower  spectral  matrices  are  written 
on  magnetic  tape. 

The  second  main  program  is  used  to  analyze  the  array  noise 
data  by  interrogating  the  crosspower  spectral  matrices.  The  analysis  options 
provided  by  this  program  are:  site  power  spectra*  component  average  power 
spectra,  multichannel  coherencies,  and  conventional  or  high- resolution  fre- 
quency-wavenumber spectra.  The  two  options  used  in  this  evaluation  are  the 
multichannel  coherence  and  high- resolution  frequency-wavenumber  spectra. 


The  multichannel  coherence  is  calculated  as  follows.  The 
elements  of  the  crosspower  spectral  matrix  corresponding  to  the  reference 
site  crosspowers  with  each  site  i to  be  used  in  the  coherence  calculation  are 
placed  in  an  array  PHI  dimensioned  (2,  40).  The  real  portion  of  the  cross- 
power is  in  row  1 and  the  imaginary  portion  is  in  row  2.  In  the  case  where 
coherence  is  to  be  calculated  between  the  reference  site  and  some  other 
data  site,  the  reference  site  autopower  is  placed  in  a scalar  PHX11  and 
the  data  site  autopower  in  scalar  PHI22.  The  coherence  squared  is  then 
calculated  as: 


COH2 


PHiq,  l)2  + PHI(2, 1 )2 
PHI11  * PHI22 


For  the  more  general  case  where  a multichannel  coherence 
is  calculated,  the  remaining  elements  of  the  crosspower  spectral  matrix  are 
placed  in  a matrix  C.  C is  factored  into  a triaugularlzed  matrix  S such  that 

C h (SH)S 

where  H indicates  the  conjugate  transpose  of  S . From  this,  an  optimum 
filter  set  FIL  is  designed  such  that 

(SH)(S)(FIL)  * PHI. 
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The  multichannel  coherence  squared  is  then  calculated  as 


coh2  = ismmM.  . 

PHIU 

Coherence  is  then  a measure  of  the  similarity  if  functions. 

For  this  evaluation  of  the  Iranian  Long-Period  Array  noise  field  characteristics, 
multichannel  coherencies  (measured  as  coherence  squared)  were  computed 
from  ten  noise  samples  using  sites  2,  4,  5,  6,  and  7 to  predict  site  1.  The 
results  were  averaged  together  to  eliminate  minor  day-to-day  variations. 

The  average  coherence  squared  values  are  plotted  in  Figures  III- 8 to  HI- 10. 

The  vertical  dashed  lines  in  each  plot  represent  the  bandpass  filter  limits 
(0.  023  - 0.  059  Ha)  used  in  RMS  noise  computations  and  signal  processing 
and  analysis.  In  each  figure,  the  coherence  squared  has  a peak  at  approxi- 
mately the  same  frequency  (0.  06  - 0.  07  Hz)  as  the  microseismic  peak  of  the 
noise  RMS  amplitude  spectra.  Robinson  (Robinson,  1967)  notes  that  micro- 
seismic  noise  is  correlated  to  various  degrees  both  in  time  and  in  space. 

This  would  explain  the  presence  of  this  peak  in  the  coherence  squared  plots. 

Within  the  signal  bandpass  filter  limits  of  0.  023  - 0.  059  Hz, 
a second  peak  appears  on  the  plots  of  coherence  squared  at  approximately 
0.  035  Hz.  The  peak  on  the  plot  of  the  vertical  component  coherence  squared 
is  lower  than  the  peaks  on  either  of  the  coherence  squared  plots  for  the  hori- 
zontal components.  The  level  of  coherence  in  the  signal  paesband  is  high 
enough  that  multi-channel  filter  processing  may  be  effective.  In  general, 
the  shape  and  coherence  level?,  shown  in  these  figures  are  quite  similar  to 
those  determined  from  the  inner-ring  sites  of  the  Norwegian  Seismic  Array 
(Laun,  Shen,  and  Swindell,  1973), 

The  frequency-wavenumber  spectra  are  calculated  in  the 
following  manner.  In  order  to  beamsteer  an  array  so  that  it  enhances  plane 
waves  from  a particular  direction,  time  delays  are  applied  to  the  data  to 
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(coherence) 


The  multichannel  coherence  squared  is  then  calculated  as 


COH 


2 _ (PHIHFIL) 
PHI  11 


Coherence  is  then  a measure  of  the  similarity  of  functions. 

For  this  evaluation  of  the  Iranian  Long- Period  Array  noise  field  characteristics, 
multichannel  coherencies  (measured  as  coherence  squared)  were  computed 
from  ten  noise  samples  using  sites  2,  4,  5,  6,  and  7 to  predict  site  1.  The 
results  were  averaged  together  to  eliminate  minor  day-to-day  variations. 

The  average  coherence  squared  values  are  plotted  in  Figures  III- 8 to  III- 10. 

The  vertical  dashed  lines  in  each  plot  represent  the  bandpass  filter  limits 
(0.  023  - 0.  059  Hz)  used  in  RMS  noise  computations  and  signal  processing 
and  analysis.  In  each  figure,  the  coherence  squared  has  a peak  at  approxi- 
mately the  same  frequency  (0.  06  - 0.  07  Hz)  as  the  microseismic  peak  of  the 
noise  RMS  amplitude  spectra.  Robinson  (Robinson,  1967)  notes  that  micro- 
seismic  noise  is  correlated  to  various  degrees  both  in  time  and  in  space. 

This  would  explain  the  presence  of  this  peak  in  the  coherence  squared  plots. 


Within  the  signal  bandpass  filter  limits  of  0,  023  - 0.  059  Hz, 
a second  peak  appears  on  the  plots  of  coherence  squared  at  approximately 
0.  035  Hz.  The  peak  on  the  plot  of  the  vertical  component  coherence  squared 
is  lower  than  the  peaks  on  either  of  the  coherence  squared  plots  for  the  hori- 
zontal components.  The  level  of  coherence  in  the  signal  passband  is  high 
enough  that  multi-channel  filter  processing  may  be  effective.  In  general, 
the  shape  and  coherence  levels  shown  in  these  figures  are  quite  similar  to 
those  determined  from  the  inner- ring  sites  of  the  Norwegian  Seismic  Array 
(Laun,  Shen,  and  Swindell,  1973). 


The  frequency- wavenumber  spectra  are  calculated  in  the 
following  manner.  In  order  to  beamateer  an  array  so  that  it  enhances  plane 
waves  from  a particular  direction,  time  delays  are  applied  to  the  data  to 
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FIGURE  III -8 

VERTICAL  COMPONENT  AVERAGE  MULTICHANNEL  COHERENCE  OF 

10  ILPA  NOISE  SAMPLES 
SITE  1 PREDICTED  FROM  SITES  2,  4,  5,  6,  7 
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FIGURE  HI- 9 

NORTH  COMPONENT  AVERAGE  MULTICHANNEL 
COHERENCE  OF  10  ILPA  NOISE  SAMPLES 
SITE  1 PREDICTED  FROM  SITES  2(  4,  5,  6,  7 


111- 21 


(coherence) 


1.0, 


0.  9 


Frequency  (Has) 


0—1 > - ..J t. — J I 

60  40  20  10  8 6 5 

Period  (sec.  ) 


FIGURE  IE- 10 

EAST  COMPONENT  AVERAGE  MULTICHANNEL 
COHERENCE  OF  10  ILPA  NOISE  SAMPLES 
SITE  1 PREDICTED  FROM  SITES  2,  4,  5,  6,  7 
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time-align  the  arrival  of  wavefront*  associated  with  that  direction.  The 
application  of  a time  delay  7^  to  the  i**1  channel  1*  equivalent  to  applying  a 
convolution  filter  ^(r)  B ^(T-r^)  who*e  Fourier  transform  it  exp(12ff£  •!?'), 
where!?  is  the  wavenumber  space  vector  and  X^is  the  coordinate  vector  for 
■ite  i.  The  beamsteer  filter  set  for  each  wavenumber  "E  is  specified  by  the 


vector 


ex  p^nTseS^) 

exp(i2»r£*5?  ) 

e 2 


jj»xp(i2?rk*XN) 

where  N is  the  number  of  channels. 


The  power  output  of  the  filter  set  is  then 
CS(  k ) * VT0V 


HRS(  k ) 


VT*'1  V 


where  CS  * conventional  frequency-wavenumber  spectra 

HRS  * high-resolution  frequency-wavenumber  spectra 
♦ = crosspower  spectral  matrix. 

By  incrementing  the  phase  vector  V by  appropriate  discrete 
values,  the  wavenumber  spectra  for  each  corresponding  point  in  Tc  space  is 
computed.  The  spectra  are  next  converted  to  dB  units  and  for  purposes  of 
plotting  the  maximum  value  is  assigned  the  symbol  +A.  Each  power  level 
below  this  maximum  is  assigned  a symbol  from  A to  Z,  the  dB  decrement 
being  user  specified. 
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Figure  HI- 11  shows  a high- resolution  frequency- wavenumber 
spectrum  computed  at  a frequency  of  0.  03516  Hz  using  the  vertical  component 
data  of  event  1295.  A signal  was  used  for  this  figure  to  show  how  well  the 
method  can  determine  arrival  azimuth  and  velocity.  The  azimuth  of  this 
event  as  computed  from  the  epicenter  coordinates  is  38.  2°,  From  the  fre- 
quency-wavenumber spectrum  of  Figure  111- 11#  the  azimuth  is  40.  9°. 

The  phase  velocity  measured  from  this  figure  is  3.  6 km/ecc. , which  appears 
to  be  a good  Rayleigh  wave  phase  velocity. 

In  order  to  investigate  the  question  of  propagating  noise, 
high  resolution  frequency-wavenumber  spectra  were  computed  using  each 
component  of  91  noise  samples.  These  spectra  were  computed  at  three 
frequencies:  at  0.  05859  Hz  for  the  microseismic  peak,  at  0.  04297  Hz  for 
a representative  signal-gate  frequency,  and  at  0.  01172  Hz  for  a point  on  the 
other  side  of  the  signal  gate.  The  arrival  azimuth  and  phase  velocity  of  the 
peak  value  of  each  frequency-wavenumber  spectrum  were  then  measured. 

Table  XIX-4  presents  the  number  of  occurrences  by  azimuth  of  the  peak  power 
of  these  frequency-wavenumber  spectra.  These  values  are  tabulated  without 
regard  for  the  phase  velocity  associated  with  the  measured  azimuth. 

The  values  in  this  table  do  not  show  any  particular  pattern  as 
regards  arrival  azimuth.  This  would  imply  that  there  is  no  dominant  source 
of  propagating  noise. 

The  next  step  in  considering  the  question  of  propagating  noise  is 
to  take  into  account  the  measured  phase  velocities.  Propagating  noise  with 
phase  velocities  well  below  or  above  the  beamforming  velocities  (4.  0 km/sec. 
for  Love,  3,  5 km/sec.  for  Rayleigh)  will  be  suppressed  by  the  beamforming 
process.  Therefore,  those  peaks  in  the  frequency-wavenumber  spectra 
which  showed  phase  velocities  below  3,  2 km/sec.  or  above  4.5  km/sec. 
were  removed  from  the  results.  The  remainder  were  grouped  by  arrival 
azimuth, 
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FREQUENCY- WAVENUMBER  SPECTRUM  FOR  VERTICAL  COMPONENT  OF  EVENT  1295 

COMPUTED  AT  FREQUENCY  = 0.  03516  Hz 


The  results  are  shown  in  Figure  III -12  lor  the  frequency- 
wavenumber  spectra  computed  at  0.  01172  Hz,  in  Figure  III- 13  for  the  fre- 
quency-wavenumber spectra  computed  at  0.  04297  Hz,  and  in  Figure  HI-14 
for  the  frequency-wavenumber  spectra  computed  at  0.  05859  Hz,  The  darkened 
area  in  each  azimuth  bin  (0°-45°,  45°-90°,  and  so  on)  represents  the  number 
of  times  propagating  noise  was  observed  arriving  at  the  array  in  that  azimuth 
range  at  signal  velocities.  Bearing  in  mind  the  above  description  of  how 
these  figures  were  created,  the  following  points  should  be  noted  from  the 
figures: 

e Figure  III- 12.  There  appears  to  be  very  little  0.  01172  Hz 

propagating  noise  arriving  at  the  array  at  signal  velocities. 

What  little  there  is  shows  no  dominant  range  of  source  azimuths. 

• Figure  III- 13.  Approximately  half  of  the  frequency-wavenumber 

spectra  computed  at  0.  04297  Hz  had  peaks  with  phase  velocities 
in  the  signal  velocity  range,  indicating  that  propagating  noise 
may  well  form  a significant  part  of  the  noise  field  at  signal 
frequencies.  On  the  vertical  and  east  components,  the  majority 
of  propagating  noise  arrivals  had  arrival  azimuths  between  90° 
and  225°,  indicating  noise  sources  lying  to  the  south  of  the 
array,  away  from  the  general  area  of  interest.  On  the  north 
component,  the  majority  of  propagating  noise  arrivals  had 

O 0 

arrival  azimuths  between  315  and  45  , 

e Figure  HI- 14.  Approximately  half  of  the  frequency-wavenumber 

spectra  computed  at  0.  05859  Hz  had  peaks  with  phase  velocities 
in  the  signal  velocity  range,  indicating  that  propagating  noise 
forms  a significant  portion  of  the  microseismic  noise  peak.  On 
all  three  components,  the  great  majority  of  the  propagating  noise 
arrivals  had  source  azimuths  to  the  south  of  the  array,  away 
from  the  general  area  of  interest. 
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SECTION  IV 
SIGNAL  ANALYSIS 


A.  DISCUSSION 

This  section  will  look  at  two  aspects  of  signal  analysis -gains 
in  signal-to-noise  ratio  due  to  beamforming  and  site-to-site  signal  similarity. 
The  first  aspect  is  intended  to  indicate  how  much  better  than  a single- site 
station  at  the  same  location  the  array  will  perform  in  terms  of  detection  cap- 
ability. To  carry  out  this  study,  aignal-to-noise  ratios  were  computed  for 
a suite  of  1 00  events  which  were  detected  on  both  the  reference  site  and  beam 
data.  The  signal-to-noise  ratio  differences  between  the  reference  site  data 

and  beam  data  provide  the  array  gain  due  to  beamforming. 

» 

The  site-to-site  signal  similarity  computations  are  intended  to 
describe  how  much  alike  a signal  recorded  at  one  site  is  to  the  same  signal 
recorded  at  other  sites.  Local  geologic  differences  from  site  to  site  and  (un- 
intended) differences  in  instrumentation  can  both  affect  the  signal  similarity. 
The  less  similar  the  signals  are  from  site  to  site,  the  poorer  the  beamforming 
process  will  work,  since  the  dissimilar  portions  of  the  signals  will  add  random- 
ly. 


B.  SIGNAL-TO-NOISE  RATIO  GAINS  DUE  TO  BEAMFORMING 

When  the  data  recorded  at  the  individual  sites  of  an  array  are 
formed  into  beams,  the  signal-to-noise  ratio  of  each  component  is  increased 
due  to  suppression  of  noise.  In  the  ideal  case,  the  noise  is  purely  random  and 
is  suppressed  by  a factor  approximately  equal  to  the  square  root  of  the  number 
of  sites  used  in  the  beamforming  process.  In  practice,  the  noise  is  composed 
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of  a random  element  and  a propagating  non- random  element.  This  propagating 
element  is  suppressed  to  a lesser  degree  than  the  random  element,  the 
amount  of  suppression  depending  on  how  far  off  the  beamforming  azimuth 
its  azimuth  lies  and  how  far  from  signal  velocities  is  its  velocity.  Also, 
the  beamforming  process  suppresses  the  signal  to  some  extent.  This 
is  dependent  on  how  accurate  the  computed  time  delays  used  to  time-align 
the  individual  traces  are  and  how  similar  the  signals  are  from  site  to  site. 

In  particular,  at  some  point  close  to  the  array,  the  plane-wave  assumption 
used  to  compute  these  time  delays  must  break  down. 

To  obtain  an  estimate  of  the  signal-to-noise  ratio  gains  to  be 
expected  from  the  beamforming  process,  a suite  of  events  was  selected  which 
were  detected  on  both  the  reference  site  and  beam  traces  and  contained  only 
noise  in  the  time  gate  immediately  preceding  the  signal  arrival  time.  The 
signal-to-noise  ratios  for  all  components  of  the  reference  site  and  beam 
traces  were  then  computed  using  the  equation: 

S/N  ,dB)  = 20,  .MO,, 

where  'zero-to-peak  amplitude1  is  the  amplitude  of  the  largest  peak  of  0 " 
signal  waveform  and  'RMS  noise'  is  measured  in  the  time  gate  immedi-i^iy 
preceding  the  signal  arrival.  The  gain  due  to  beamforming  is  then  simply  the 
difference  between  the  beam  signal-to-noise  ratio  and  the  reference  site  signal- 
to-noise  ratio. 

The  results  in  Table  IV- 1 are  grouped  on  the  basis  of  epi- 
central  distance.  The  values  in  the  column  headed  'optimum  gain1  were 
computed  from  the  following  equation: 

Optimum  Gain  (dB)  “20,  ’I'LOG^  yj number  of  sites 
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where  'average  number  of  Bites'  is  the  average  number  of  siteB  used  in  beam- 
forming. This  table  shows  that  beamforming  gains  for  all  components  of 
events  with  epicentral  distances  less  than  ten  degrees  are  very  low.  One  pos- 
sible explanation  of  this  is  that  the  plane-wave  assumption  used  in  beamform- 
ing to  compute  time  delays  fails  for  eventB  with  epicenters  less  than  ten  de- 
grees from  the  array.  A second  possible  explanation  is  that  the  fixed  veloc- 
ities used  in  beamforming  (3.  5 km/ sec  for  Rayleigh,  and  4.  0 km/sec  for 
Love)  are  not  appropriate  for  beamforming  close  events. 

The  mean  gains  for  the  other  ranges  of  epicentral  distances 
remain  fairly  constant.  This  Implies  that  the  plane- wave  assumption  holds 
for  events  with  epicentral  distances  greater  than  ten  degrees. 

Comparing  the  mean  gains  in  Table  IV- 1 with  the  corresponding 
optimum  gains,  the  data  show  that  in  general  the  mean  gains  for  the  horizontal 
components  are  lower  than  the  optimum  gains.  This  implies  that  some  of  the 
noise  is  propagating,  since,  as  was  described  earlier,  propagating  noise 
is  suppressed  by  beamforming  to  a lesser  degree  than  is  random  noise. 

An  interesting  feature  of  the  data  in  Table  IV- 1 is  that  the 
radial  component  gains  are  lower  than  the  vertical  component  gains.  In  Table 
IV-2  the  signal-to-noise  ratio  gains  are  separated  into  the  signal-to-noise  gain 
due  to  RMS  noise  suppression  and  the  signal-to-noise  ratio  loss  due  to  peak 
signal  suppression.  From  the  data  in  this  table  it  appears  that  the  difference 
in  signal-to-noise  ratio  gain  between  the  vertical  and  radial  components  is  due 
roughly  equally  to  both  lower  RMS  noise  suppression  and  greater  signal 
suppression  on  the  radial  component.  The  lower  RMS  noise  suppression  on 
the  radial  component  in  conjunction  with  lower  RMS  noise  suppression  on 
the  transverse  component  relative  to  the  vertical  component  implies  that  there 
is  more  propagating  noise  on  the  horizontal  components  than  on  the  vertical 
components  and  that  this  noise  is  propagating  as  both  Love  and  Rayleigh  waves. 
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Epicentral 

Distance 


0°-10°  6.5  5.0  5. 


10°-20°  8.8  5.8 


TABLE  IV- 2 

NOISE  AND  SIGNAL  SUPPRESSION  IN  dB 
DUE  TO  BEAMFORMING 


RMS  NoLae  Suppression  Peak  Signal  Suppression 


C.  SITE-TO-S1TE  SIGNAL  SIMILARITY 

The  program  used  to  measure  site-to-site  signal  similarity 
requires  as  its  input  the  edited  time-domain  data  in  the  recorded  vertical, 
north,  east  configuration.  The  program  computes  beamsteer  time  delays 
for  each  propagation  mode  by  the  equation: 

(X  -X.)  sin  0 + (Y  -Y.)  cos  $ 

n ss  — E — ± E i 

1 V 

where:  is  the  time  delay  for  site  i 

X and  Y are  the  X and  Y coordinates  of  the  reference  site 
r r 

Xi  and  Yj  are  the  X and  Y coordinates  of  the  data  site  i 

0 is  the  beamsteer  azimuth,  and 

Y is  the  velocity  of  the  propagation  mode. 

The  time  delay,  or  lag,  is  considered  to  be  negative  if  the 
signal  arrives  at  data  site  i before  it  arrives  at  the  reference  site  and  positive 
if  it  arrives  at  the  data  site  i after  it  arrives  at  the  reference  site. 

The  data  at  the  various  sites  next  have  their  trace  means  re- 
moved and  are  rotated  to  the  beamsteer  azimuth.  Following  this,  the  reference 
site  zero-lag  autocorrelation  values  are  computed  from: 

*rr<°>  4 W1* 

where  X^  (k)  represents  the  reference  site  time  series  of  N points  sampled  at 

dt  intervals.  Correlation  processing  then  continues  for  the  data  sites.  The 

zero-lag  autocorrelation  $^(0)  is  first  computed.  The  cross-correlation 

functions  0 At  A t)  are  next  computed  in  the  user- specified  range  for  each 
ri 

lag  T.  The  zero -lag  center  point  of  the  range  of  lags  for  a site  is  the  beamsteer 
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lag  for  that  site.  The  cross -correlation  function  for  each  site  i at  a lag  r 
is  computed  using 


*xi(TAt)  = ~ xt(k)  Xt(k+r4t) . 

After  the  autocorrelation  values  and  cross- correlation 
functions  are  computed  for  each  mode -component  combination  at  each  site, 
the  cross- correlation  matrix  is  searched  for  its  maximum  values  and 
corresponding  lags.  For  each  combination  the  correlation  coefficient  is 
computed  according  to 

CC  - **iU) 

ri  orr<°>*11«»]1 

■where  A is  the  lag  at  which  ^(T^t)  is  a maximum  and  CC^  is  the  correla- 
tion coefficient  for  the  reference  site  and  date  site  i.  Correlation  coefficient 
means  are  computed  after  all  sites  are  processed. 


s- 


Site-to-site  signal  similarity  was  investigated  for  each 
component  by  generating  correlation  coefficients  between  the  reference  site 
and  the  remaining  six  sites  for  a suite  of  large  events.  The  results  of  com- 
puting correlation  coefficients  for  22  events  ranging  in  bodywave  magnitude 
from  4.  7 to  6.  3 and  one  typical  noise  sample  are  shown  in  Tables  IV-3  to 
1V-5  for  the  vertical,  transverse,  and  radial  components  of  motion,  respect- 
ively. Included  in  each  table  are  the  correlation  coefficient  for  each  site-event, 
the  signal-to-noise  ratio  of  the  reference  site  data  for  each  event,  the  average 
correlation  coefficient  for  each  event,  and  the  correlation  coefficient  for  each 
site  averaged  over  all  events.  A dash  indicates  that  data  for  that  site-event 
were  either  not  available  or  were  not  suitable  for  processing.  No  signal-to- 
noise  ratio  is  listed  for  events  whose  noise  gates  contained  other  signals  or 
spikes.  Since  no  plane-wave  anomalies  were  observed,  cross -cor relation 
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TABLE  IV -3 

CORRELATION  COEFFICIENTS  FOR  VERTICAL  COMPONENT 


Event 

Number 

Reference 

Site 

S/N  (dB) 

2 

3 

H 

6 

7 i 

Event  Average 
Correlation 
Coefficient 

0766 

36.  2 

0.80 

m m m m 

0.83 

0.  65 

0.  80 

0.95 

0.81 

0790 

32.  8 

• . . 

.... 

0.  76 

0.  83 

0.  84 

0.85 

0.82 

0824 

19.  6 

0.  54 

--  -- 

0.  56 

0.  65 

0.77 

0.20 

0.  53 

0885 

— 

0.99 

— 

— 

0.  95 

0.  95 

0.  95 

0.96 

0886 

60.2 

0.84 

— 

0.  82 

0.82 

0.  93 

0.  78 

0.84 

0890 

36.  7 

0.89 

m m „ m 

0.  84 

0.87 

0.  94 

0.81 

0.87 

0900 

38.  0 

0.84 

— 

0.84 

0.  80 

0.  95 

0.78 

0.84 

1265 

— 

0.  51 

— 

0.75 

0.  71 

0.33 

0.57 

1295 

31.4 

0.  93 

0.  90 

0.93 

0.  91 

0.89 

0.91 

, 1296 

74.4 

0.81 

— 

0.  58 

0.  74 

0.86 

0.  92 

0.  78 

45  1 

0 

0.  84 

0.  81 

0 83 

JL  Jwl 

1331 

55.  1 

0.80 

..... 

.... 

.... 

0.88 

0.70 

0.79 

1395 

34.  9 

0.85 

.... 

0.  77 

0.  90 

0.91 

0.86 

1406 

39.  5 

0.85 

0.  92 

.... 

0.  94 

0.  95 

0.82 

0.90 

1413 

50.4 

0.  91 

0.  96 

— 

0.  92 

— 

0.87 

0.  91 

1524 

56.  2 

0.86 

0.  90 

0.  78 

.... 

0.91 

.... 

0.86 

1544 

33.  9 

0.  92 

0.  95 

.... 

.... 

— 

0.82 

0.90 

1555 

35.4 

0.84 

0.85 

.... 

0.  78 

0.89 

0.64 

0.80 

1574 

30.5 

0.76 

0.87 

— 

.... 

0.82 

0.  62 

0.77 

1621 

54.4 

— 

0.88 

0.  83 

.... 

0.72 

0.81 

1625 

33. 1 

. m m m 

MMMM 

0.  87 

0.  84 

0.  91 

0.78 

0.85 

1628 

31.  8 

— 

— 

0.  82 

0.  94 

0.75 

— 

0.  83 

Day  170 

(1976) 

— 

0.26 

— 

0.  17 

-0.  15 

-0. 14 

-0.19 

-0.08 

Noise  Sample 

Site  Average 

(Excluding 

0.  86 

0.  90 

0.  81 

0.  84 

0.88 

0.81 

.... 

Events  824 

and  1265) 
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TABLE  IV-4 

CORRELATION  COEFFICIENTS  FOR  TRANSVERSE  COMPONENT 


Event 

Number 

Reference 

Site 

S/N  <dB) 

2 

3 

4 

5 

6 

7 

Event  Average 
Correlation 
Coefficient 

0766 

41.3 

0.81 

m — m 

0.87 

0.89 

0.  96 

0.  92 

0.89 

0790 

38.7 

mmmm 

0.75 

0.  66 

0.  94 

0.  95 

0.88 

0824 

22.0 

0.  71 

mmmm 

0.  54 

0.65 

0.  78 

0.  62 

0.  66 

0885 

m m m m 

0.  82 

M M • - 

.... 

0.  86 

0.  95 

0.79 

0.86 

0886 

47.  1 

0.81 

0.  81 

0.80 

0.  97 

0.  82 

0.B4 

0890 

36.7 

0.89 

mmmm 

0.  84 

0.87 

0.  96 

0.88 

0.89 

0900 

41.5 

0.82 

mmmm 

0.  75 

0.76 

0.  94 

0.  82 

0.82 

1265 

m m m w 

0.57 

mmmm 

0.  62 

0.55 

.... 

0.  60 

0.  58 

1295 

37.0 

0.  92 

mmmm 

0.76 

0.  92 

0.84 

0.  86 

0.86 

1296 

68.3 

0.85 

— 

0.81 

0.  84 

0.89 

0.85 

0.  85 

1321 

50.4 

0.94 

. . . . 

««••■■■ 

mmmm 

0.  96 

0.  94 

0,  95 

1331 

45.6 

0.85 

---- 

mmmm 

.... 

0.  93 

0.  83 

0.  87 

1395 

39.0 

0.78 

.... 

mmmm 

0.88 

0.  93 

0.  88 

0.87 

1406 

36.  3 

0.  72 

0.89 

.... 

0.  68 

0.89 

0.  66 

0.77 

1413 

47.  0 

0.  93 

0,  96 

..... 

0.91 

.... 

0.89 

0.  92 

1524 

57.  5 

0.  85 

0.  89 

0.  86 

.... 

0.  91 

mmmm 

0.  88 

1544 

45. 1 

0.  97 

0.  86 

— 

— 

. — 

0.95 

0.  92 

1555 

30.4 

0.87 

0.  87 

— 

0.  82 

0.  92 

0.  96 

0.  97 

1574 

49.9 

0.  97 

0.  98 

wm  m m m 

.... 

0.  98 

0.80 

0.86 

1621 

60.3 

m>  ■§  m * 

0.97 

0.87 

.... 

.... 

0.  92 

0.  92 

1625 

28.8 

m m m — 

m m m m 

0.  65 

0.  68 

0.  93 

0.72 

0.  74 

1628 

34.2 

— 

mmmm 

0.  75 

0.  95 

0.58 

.... 

0.76 

Day  170 

(1976) 

— 

0. 10 

mmmm 

0.  06 

-0.14 

0.21 

0.21 

0.  09 

Noise  Simple 

Site  Average 

(Excluding 

0.86 

0.  92 

0.7S 

0.84 

0.  91 

0.86 

•*•*“»* 

Eventi  824 

and  1265) 

r.  il.j- 
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lag  data  are  omitted.  The  data  in  these  tables  are  interpreted  from  the  stand- 
point that  a correlation  coefficient  of  1 computed  using  data  recorded  at  the 
reference  site  and  data  site  i would  indicate  that  the  waveforms  at  these  sites 
are  identical,  a correlation  coefficient  of  zero  would  indicate  that  the  wave- 
forms at  these  sites  are  wholly  unlike  each  other,  and  a correlation  coefficient 
of  -1  would  indicate  that  the  waveforms  at  these  sites  differ  by  a phase  shift 
of  180°. 

The  following  points  should  be  noted  from  this  data: 

d 

• Events  with  low  correlation  coefficients  were  reprocessed 
using  site  2 as  reference  to  determine  if  site  1 was  producing 
anamolous  values.  The  results  were  not  significantly  different.  1 

• The  low  correlation  coefficients  for  event  824  appear  to  be  a 

result  of  low  signal-to-noise  ratio.  Note  that  this  event  has  1 

the  lowest  signal-to-noise  ratio  of  the  suite  of  events. 

• Event  1265  also  displayed  low  correlation  coefficients.  Plots 
of  the  raw  data  of  this  event  showed  the  data  to  be  corrupted 
with  high  frequency  noise.  Examination  of  the  field  tape 
logs  for  this  time  period  revealed  that  transmission  errors 
occurred  during  the  recording  of  this  event. 

• The  transverse  component  yielded  larger  correlation  co-  |< 

efficients  than  either  the  vertical  or  radial  components,  being 

greater  than  0.  8 for  all  large  events. 

• The  radial  component  yielded  on  average  the  lowest  correlation 
coefficients.  This  is  perhaps  reflected  in  the  results  of  peak 
signal  suppression  due  to  beamforming  the  radial  component 
as  shown  in  Table  IV- 2. 
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Overall,  sites  3 and  6 yield  better  correlation  coefficients 
with  the  reference  site  than  do  the  other  sites.  Site  4 on  the 
average,  yields  the  lowest  correlation  coefficients. 

As  expected,  the  noise  sample  yielded  very  low  correlation 
coefficients. 

The  peak  signal  suppressions  noted  in  Table  IV- 2 appear 
to  be  reflected  by  the  less  than  perfect  site-to-Bite  sim- 
ilarities. 
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SECTION  V 

ILPA  DETECTION  CAPABILITY 


A.  DISCUSSION 

In  past  evaluation  tasks,  detection  statistics  were  derived  from 
an  event  population  for  every  member  of  which  a clear  detection /non-detection 
decision  could  be  made;  i.  e. , the  analyst  could  state  that  he  either  saw  the 
Kought-for  signal  or  he  saw  seismic  noise.  Unfortunately,  the  world  does  not 
always  present  the  analyst  with  such  a clear-cut  case.  Mixed  events,  system 
failures  resulting  in  no  data  being  recorded,  and  malfunctions  all  tend  to 
obscure  the  detection  capability  picture. 

The  term  'mixed  event*  refers  to  the  case  where  the  sought- 
for  signal  is  obscured  or  completely  masked  by  a second  signal.  This  can 
happen  either  when  the  two  signals  arrive  at  the  station  at  essentially  the  same 
time  or  when  a larger  signal  arrives  before  the  signal  under  analysis,  burying 
this  signal  in  its  coda.  The  term  'system  failure*  refers  to  the  total  shutting 
down  of  the  station  so  that  no  data  is  recorded.  The  term  'malfunction'  refers 
to  any  partial  failure  of  the  system,  from  sensor  unit  to  reception  of  data  at  the 
Seismic  Data  Analysis  Center,  which  causes  degradation  of  the  seismic  data. 

The  problem  of  mixed  events  1b  often  difficult  for  the  analyst 
to  resolve  and  is  probably  the  major  source  of  false  alarms.  (The  term 
'false  alarm'  in  this  context  means  declaring  a detection  when  in  fact  the 
observed  signal  is  from  an  event  other  than  that  under  analysis. ) When  a 
signal  is  observed  in  the  time  gate  of  the  event  under  analysis,  the  analyst 
first  checks  the  waveforms  on  the  three  components  of  motion  to  see  that  their 
inter-relationships  are  correct.  If  still  in  doubt,  the  analyst  checks  available 
event  lists  to  see  whether  any  other  reported  event  could  have  arrived  in  the 
signal  gate.  In  general,  the  analyst  declares  a detection  if  a dispersed  signal 
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is  observed  having  the  correct  inter-relationships  between  the  Love  and 
Rayleigh  waves  and  if  no  other  event  has  been  reported  which  could  be  mis- 
taken for  the  event  under  analysis. 


I 

I 
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Since  mixed  events  and  (secondarily)  events  for  which  no  data 
were  recorded  or  which  contained  malfunctions  are  a fairly  common  problem, 
the  Iranian  Long- Period  Array  detection  capability  estimates  are  calculated 
in  two  ways.  The  first  of  these  is  termed  the  absolute  detection  capability 
estimate.  When  computing  this  estimate,  all  mixed  events,  events  for  which 
no  data  were  recorded,  and  events  containing  malfunctions  are  counted  as 
non -detections  when  forming  the  detection  statistics.  This  approach  gives 
a real-world  detection  capability  estimate. 

The  second  of  these  estimates  is  termed  the  conditional  detec- 
tion capability  estimate.  When  computing  this  estimate,  all  mixed  events, 
events  for  which  no  data  were  recorded,  and  events  containing  malfunctions  are 
rejected  from  the  detection  statistics.  This  approach  gives  an  ideal  detection 
capability  estimate,  The  value  of  this  ideal  estimate  is  that  it  shows  the  de- 
tection capability  improvement  possible  if  the  reliability  of  the  Instrumentation 
can  be  improved  and  if  methods  of  separating  mixed  events  can  be  found. 

The  number  used  to  represent  detection  capability  is  the  50 
percent  detection  threshold,  denoted  by  The  50  percent  detection 

threshold  is  the  bodywave  magnitude  for  which  the  probability  of  detection  is 
0.  5.  It  is  computed  by  fitting  the  Gaussian  probability  function  to  the  detec- 
tion statistics  by  a maximum  likelihood  method  (Ringdal,  1974).  Hereafter, 
this  will  be  referred  to  as  the  maximum  likelihood  curve. 

It  must  be  kept  in  mind  that,  since  the  data  base  was  derived 
from  the  NORSAR  and  NEIS  event  bulletins,  all  detection  capability  estimates 
are  in  terms  of  a combined  'NORSAR  - NElS'm,  unit.  This  fact  is  Important, 
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since  NORSAR  units  are  not  the  same  as  those  from  the  National  Earth- 
quake Information  Service  (NEIS)  event  lists. 

Depending  on  the  quantity  of  data  available,  one  of  three  levels 
of  confidence  is  placed  on  each  50  percent  detection  threshold.  If  a 50  per- 
cent detection  threshold  is  given  to  two  decimal  places  (e.  g. , = 4.  56), 

sufficient  detection  statistics  were  available  to  compute  a reliable  detection 
capability  estimate.  If  a 50  percent  detection  threshold  is  enclosed  in  par- 
entheses and  given  to  only  one  decimal  place,  (e.  g. , = (*•  &))  the 

detection  statistics  were  sparse  and  the  detection  capability  estimate  can 
be  considered  to  be  only  a first  approximation.  Finally,  If  a 50  percent  de- 
tection threshold  is  expressed  as  greater  than  some  value  (e.  g, , m^>5.0), 
the  detection  statistics  were  too  sparse  to  allow  fitting  of  a maximum  like- 
lihood curve.  In  most  cases,  this  is  due  to  a lack  of  detected  events.  When 
this  occurred,  the  value  given  is  the  m^  of  the  largest  non-detected  event. 

H.  ILPA  LONG-PERIOD  DETECTION  CAPABILITY  ESTIMATES 

This  subsection  examines  the  estimates  of  detection  capability 
made  from  the  accumulated  detection  statistics  for  the  reference  site  and 
beam  data.  The  criteria  used  by  the  analyst  which  determine  whether  an 
event  was  detected  are! 

• The  presence  of  dispersion  in  the  signal  gate. 

e A peak  in  the  dispersed  wavetraln  3 dB  or  more  above  any 

peak  outside  the  dispersed  wavetrain  and  inside  a time  gate 
starting  600  seconds  before  the  predicted  Love  wave  arrival 
time  and  ending  600  seconds  after  the  estimated  Rayleigh 
wave  end  time. 

e Occurrence  of  signal  onset  within  ±180  seconds  of  the  pre- 

dicted signal  onset  time. 
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• Detection  of  the  signal  on  at  least  two  of  the  three  components 

of  motion. 

These  criteria  were  used  as  a guide  to  aid  the  analyst  In 
determining  the  detection  status  of  processed  events.  The  first  criterion 
was  always  followed  for  the  events  processed  in  this  evaluation,  since  this 
is  the  primary  visual  difference  between  seismic  signals  and  noise.  The 
second  criterion  was  Occasionally  disregarded,  since  bodywaves  such  as 
SS  or  isolated  noise  pulses  would  be  at  times  visible  within  the  specified 
time  gate.  If  the  event  under  analysis  was  visible,  it  was  declared  detected 
even  if  one  of  these  was  also  visible.  The  third  criterion  was  also  occasion- 
ally not  followed.  When  a waveform  in  the  signal  gate  was  observed  to  arrive 
later  than  this  criterion  allowed,  the  event  lists  were  checked  to  ensure  that 
this  waveform  was  not  due  to  some  other  event.  If  no  other  event  could  be 
found  whose  surface  waves  would  arrive  at  the  observed  arrival  time,  the 
event  under  analysts  was  called  a detection.  The  last  criterion  was  rigidly 
followed,  since  it  was  Imposed  to  reduce  the  probability  of  erroneously  de- 
claring an  event  to  be  detected. 

The  various  detection  capability  estimates  made  for  this  eval- 
uation are  shown  in  Figures  V-l  to  V-4.  Each  'sub-figure'  of  these  figures 
consists  of  two  parts.  The  upper  part  consists  of  a histogram  showing  the 
number  of  detected  and  non-detected  events  as  a function  of  body-wave  mag- 
nitude (rr^)  for  the  particular  data  subset  under  consideration.  The  lower 
part  shows  the  detection  probability  derived  from  these  detection  statistics 
ar  a function  of  bodywave  magnitude.  The  percentage  of  events  detected  at 
each  m^  value  is  represented  by  an  asterisk.  The  maximum  likelihood  curve 
fitted  to  these  detection  percentages  is  represented  by  a solid  line.  The  90 
percent  confidence  limits  of  this  curve  are  represented  oy  dashed  lines. 

The  values  for  'MB50'  and  'MB90'  shewn  on  the  figures  represent  the  50 
and  90  percent  detection  thresholds  respectively,  as  determined  by  the 
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FIGURE  V-l 

ILPA  LONG- PERIOD  REFERENCE  SITE  AND 
BEAM  DETECTION  STATISTICS 
ALL  REGIONS  A AND  E 
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maximum  likelihood  curve.  Finally,  the  value  given  for  'SIGMA'  is  the 
standard  deviation  o£  the  Gaussian  probability  function  for  the  maximum 
likelihood  curve. 

Figure  V-l  presents  the  reference  site  and  beam  detection 
capability  estimates  using  the  detection  statistics  of  all  the  regions.  These 

i 

were  made  In  two  ways.  In  the  first,  denoted  by  'ALL  REGIONS  A',  mixed 
events,  events  for  which  no  data  were  recorded,  and  events  containing  mal- 
functions  were  counted  p.s  non-detections . This  forms  the  previously  defined 
absolute  detection  capability  estimate.  In  the  second,  denoted  by  'ALL 
REGIONS  B',  mixed  events,  events  for  which  no  data  were  recorded,  and 
events  containing  malfunctions  were  rejected  from  the  detection  statistics. 

This  forms  the  previously  defined  conditional  detection  capability  estimate. 

Figures  V-2  and  V-3  show  the  detection  capability  estimates 
for  the  beam  data  on  a regionalized  basis.  These  estimates  are  conditional 
detection  capability  estimates.  The  paucity  of  detection  statistics  makes  these 
regionalized  estimates  valid  as  first  approximations  only. 

The  detection  statistics  of  Figures  V-l  to  V-3  were  derived 
from  earthquakes  whose  event  parameters  were  taken  from  either  the  NORSAR 
or  the  NEIS  event  bulletins.  Since  the  two  event  bulletins  partially  overlap 
in  the  time  frame  each  covers,  it  was  possible  to  compare  the  m.  values  that 

D 

each  reports,  A total  of  518  Eurasian  events  were  found  to  be  reported  by  both 
event  bulletins.  By  fitting  a straight  line  to  the  NORSAB  - NEIS  m^  pairs  so 
found,  using  an  algorithm  which  treats  neither  variable  as  dependent,  the 
following  relationship  between  the  two  types  of  m,  was  derived: 

D 

NORSAR  m^  = 1.11  NElSm^  - 0.71  (variance  ■ 0.03), 

Figure  V-4  was  formed  by  converting  all  NEIS-reported  m^  values  In  the 
data  base  to  NORSAR  m^  values  using  the  above  relationship.  Like  Figure  V-l, 
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this  figure  shows  the  single- site  and  beam  detection  capability  estimates 
using  'ALL  REGIONS  A'  and  'ALL  REGIONS  B'  detection  statistics.  The 
purpose  of  this  is  to  estimate  the  effects  of  the  m^  conversion  on  the  de- 
tection capability  estimates. 


The  various  detection  capability  estimates  are  summarized  in 
Table  V-l.  Since  the  prime  purpose  of  this  section  is  to  estimate  the  array 
detection  bapability,  no  regionalized  reference  site  detection  capability  es- 
timates were  made.  The  following  points  should  be  noted  from  the  data  pre- 
sented in  Table  V-l: 

e Mixed  events,  events  for  which  no  data  were  recorded,  and 

malfunctions  raised  the  50  percent  detection  threshold  by 
0. 4 - 0,  5 m,  units. 


e 


The  beamforming  process  lowered  the  50  percent  detection 
threshold  by  0.  25  - 0.  30  m^  units. 


The  conversion  of  NEIS  m,  values  to  NORSAR  m^  values  had 
essentially  no  effect  on  the  ILPA  detection  capability. 


e The  regionalized  ILPA  detection  capability  estimates  can  be 

considered  as  first  approximations  only,  due  to  the  limited 
detection  statistics  available. 


It  should  be  noted  that,  although  the  conversion  of  NEIS 
values  to  NORSAR  values  had  essentially  no  effect  on  these  detection 
capability  estimates,  the  authors  feel  that  in  any  future  work  this  matter 
should  be  considered  again.  For  the  data  base  used  in  this  report,  the  ratio 
of  NORSAR  eventB  to  NEIS  events  was  approximately  4 to  1.  Future  in- 
creases in  the  data  base  will  raise  the  number  of  NEIS  events  in  the  data 
base  (since  the  NORSAR  event  bulletin  ceased  being  issued  30  September  1976) 
and  hence  will  Increase  the  relative  effect  of  NEIS  m^  values  on  the  detection 
capability  estimates. 
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TABLE  V-l 

ILPA  LONG-PERIOD  DETECTION  CAPABILITY 


Region 

NORSAR  & NEIS  rr^'s 

NORSAR  m 'a 

D 

Moan  Distance 
in  Degrees 

Reference 

Beam 

Reference 

Beam 

All  Regions  A 

4.71 

4.46 

4.  67 

4.42 

36.5 

Ail  Regions  B 

4.30 

3.  98 

4.  26 

3.  94 

J 

— 

(4.5) 

— 

(4.5) 

71.8 

2 

(4.5) 

— 

(4.4) 

69.6 

5 

— 

>4.  0 

— 

>4.  3 

40.7 

6 

— 

(3.8) 

— 

(3.7) 

36.  9 

7 

— 

(3.7) 

— 

(3.7) 

31.6 

9 

— 

(3.8) 

— 

(3.8) 

8.6 

10 

— 

(3.7) 

— 

(3.7) 

17.2 

11 

— 

(3.6) 

— 

(3.6) 

24.7 

12 

’ W a*  m m 

(3.8) 

- — — . 

(3.8) 

53,  9 

Note: 

1.  Detection  capability  1b  estimated  in  terms  of  the  50  percent 
detection  threshold. 

2.  All  Regions  A - Mixed,  no  data  recorded,  and  malfunctions  are 
counted  as  non- detections. 

3.  All  Regions  B - Mixed,  no  data  recorded,  and  malfunctions  are 
rejected  from  detection  statistics. 
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SECTION  VI 

EARTHQUAKE -PRESUMED  EXPLOSION  DISCRIMINATION 

■i 

j 

A.  DISCUSSION 

This  section  considers  the  question  of  discriminating  between 
earthquakes  and  presumed  nuclear  explosions  using  long-period  data.  The 
discrimination  method  used  is  the  surface-wave  magnitude  (M^)  versus  bodywave  j 

magnitude  (m^)  plot.  This  plot  of  Mb  versus  is  expected  to  function  as  j 

an  earthquake-presumed  nuclear  explosion  discriminant  since,  for  a given  ] 

m^,  an  explosion  generates  much  lower  Rayleigh  and  Love  waves  than  docs  j 

an  earthquake.  (Theoretically,  an  explosion  should  generate  no  Love  waves,  jj 

since  the  source  is  completely  compressional.  However,  some  Love  wave  j 

energy  is  radiated  from  tectonic  strain  release  (Sun,  1?77)  subsequent  to  the  1 

.1 

explosion.  ) Therefore,  a plot  of  M versus  m,  can  be  expected  to  show  a 

J D 

separation  of  the  data  points  into  an  earthquake  population  and  an  explosion 
population . 

The  data  base  for  this  discrimination  study  is  comprised  of 
all  events  which  were  visually  detected  on  bandpass- filtered  (0.  023  - 0.  059  Hz 
passband)  plots. 


B.  COMPUTATION  OF  SURFACE  WAVE  MAGNITUDES 


In  earlier  evaluation  tasks,  the  signal  amplitude  and  period 
values  used  in  computing  surface  wave  magnitudes  (M^)  were  measured  man- 
ually on  the  filtered  signal  plots  at  various  periods.  The  computation  of 
M^  then  uy  sd  the  aquation: 


M 


s 


LOG!0 


_ a*sf 
T*Q*G 


« LOGj  I-  1.1? 
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where: 


A = peak-to-peak  amplitude  measured  in  inches  on  the  plot. 

SF  = plot  scale  factor  in  computer  counts  per  inch, 

T = period  in  seconds  of  the  measured  amplitude. 

Q » quantization  factor  ( 20.  951  computer  counts  per  milli- 
micron). 

G **■  Instrument  response  correction  factor,  and 
A = epicentral  distance  in  degrees. 

This  approach  has  several  disadvantages.  First,  the  analyst 
is  apt  to  have  difficulties  finding  the  largest  peak  at  each  desired  period  - it 
is  sometimes  necessary  to  measure  the  periods  of  a number  of  waveforms 
before  finding  the  desired  period.  Next,  having  found  the  waveform  for  which 
Mg  is  to  be  measured,  the  analyst  may  make  a measurement  error.  Finally, 
errors  are  apt  to  occur  either  in  transcribing  the  measured  values  or  in 
calculating  the  surface-wave  magnitude  from  these  values. 

In  order  to  avoid  these  problems  and  make  the  process  of  ob- 
taining M values  less  tiresome,  a program  was  written  to  perform  the 
s 

measurements  and  calculations  automatically.  This  program  operates  by  first 
finding  the  times  of  all  zero  crossings  in  a given  time  gate  and  the  maximum 
absolute  amplitude  in  millimicrons  between  each  pair  of  adjacent  zero  cross- 
ings. The  waveform  period  is  then  simply: 

T » 2* [Time  of  Zero,,,  - Time  of  Zero.  ] . 

i+1  i ■* 

The  program  reads  the  epicentral  distance  from  the  event  header.  The 
quantization  factor  and  instrument  response  corrections  are  built  into  the 
program.  (The  instrument  response  corrections  are  derived  from  the  In- 
strument response  curve  of  Figure  VI- 1. ) 
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The  program  uses  the  above  information  to  compute  a surface- 

wave  magnitude  for  each  amplitude -period  pair.  These  M values  with  their 

8 

associated  periods  and  times  of  measurement  are  printed  out.  Finally  the 

program  picks  the  largest  M value  at  each  specified  period  of  interest  and 

s 

prints  it  out.  The  analyst  now  checks  these  largest  Mg  values  against  the  plot 
of  the  data.  If  he  does  not  like  an  M value  at  a particular  period  (for  example, 
if  the  M value  appears  to  be  associated  with  a noise  pulse  in  the  signal  gate) 
he  can  select  another  value  from  the  M list  generated  by  the  program. 

By  using  this  program,  measurement  errors  have  been  elim- 
inated and  transcription  errors  greatly  reduced  (since  the  analyst  only 
writes  down  the  final  M value  and  none  of  the  intermediate  computational 

values).  Also,  the  time  required  to  arrive  at  the  M values  for  an  event 

’ s 

has  been  reduced. 

A further  benefit  to  the  evaluation  task  was  found  once  thiB 
program  was  put  into  use.  Since  measurements  were  no  longer  being  made 
directly  on  the  plot,  it  was  possible  to  reduce  the  length  of  the  plots  by  one- 
half.  When  measurements  were  made  on  the  plot,  a horizontal  scale  of  100 
seconds  per  inch  was  considered  necessary  to  minimize  measurement  errors. 
Plots  are  now  made  with  a horizontal,  scale  of  200  seconds  per  inch.  This 
reduces  the  amount  of  plot  paper  and  plot  time  by  one-half  and  produces 
more  manageable  plots. 

During  the  preliminary  Iranian  Long- Period  Array  evaluation, 
was  measured  at  periods  of  20,  30,  and  40  seconds.  However,  measure- 
able  40  second  energy  was  not  often  found  - only  half  as  many  40  second  M 

B 

values  were  measured  as  20  second  or  30  second  M values.  It  was  also 

s 

noted  that  events  close  to  the  array  have  the  majority  of  their  energy  con- 
centrated at  or  near  25  seconds  period.  For  these  reasons,  it  has  been  de- 
cided to  drop  the  measurement  of  40  second  M and  add  the  measurement  of 
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25  second  M , 
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C.  DISCRIMINATION  RESULTS 

Table  VI- 1 lists  those  events  of  the  data  base  which  may 
be  termed  presumed  nuclear  explosions.  They  were  selected  on  the  basis 
of  their  epicentral  locations  and  bodywave  magnitudes.  These  events  were 
selected  because  their  epicenters  are  at  or  very  close  to  the  eastern  Kazakh 
test  site  and  their  bodywave  magnitudes  are  larger  than  most  earthquakes 
from  that  area. 

Figures  VI- 2 to  VI- 14  show  the  M^-m^  plots  for  the  Iranian 
Long- Period  Array  data.  The  values  are  plotted  for  the  array  with  three 
plots  encompassing  data  from  all  regions  showing  surface  wave  magnitude 
measured  at  20,  25,  and  30  seconds,  followed  by  individual  plots  for  25-second 
surface  wave  magnitudes  from  each  region.  The  symbols  used  in  these  plots 
are: 

O - earthquake  with  depth  less  than  60  km  or  depth  unknown. 

4-  - earthquake  with  depth  greater  than  60  km. 

* - presumed  nuclear  explosion  from  Region  8. 

Such  depth  information  as  was  available  came  from  the  National 
Earthquake  Information  Service  bulletins.  The  straight  line  in  each  plot  re- 
presents the  equation  of  the  M^-rn^  relationship  for  that  data  set.  This  was 
computed  using  the  data  points  for  earthquakes  with  depth  less  than  60  km 
or  depth  unknown.  The  relationship  is  computed  with  a linear  fitting  al- 
gorithm which  treats  neither  variable  as  dependent.  The  earthquakes  with 
depths  known  to  be  greater  than  60  km  were  excluded  from  the  fitting  proce- 
dures to  avoid  biassing  the  fit.  (Deeper  events  tend  to  generate  relatively 
lower  M(  values.  ) The  slopes  and  intercepts  of  the  lines  fitted  to  the  20,  25, 
and  30  second  data  are  given  in  Table  VI-2.  The  lines  fitted  to  the  data  of 
the  individual  regions  are  based  for  the  most  part  on  rather  sparse  populations 
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REGI< 


SLOPE  and 


Period 

Component 

a 

20 

Vertical 

1.51 

Transverse 

1.51 

25 

Vertical 

1.61 

Transverse 

1.49 

30 

Vertical 

1.58 

Transverse 

1.56 

where  Ma=am^+b 
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ABLE  VI- 2 
[s_mb 

INTERCEPT  VALUES 


Center  of  Mass 


b 

.V 

mb 

M. 

8 

n 

-3.  34 

0.20 

4.73 

3.83 

138 

-3.09 

0.21 

4.  64 

3.89 

152 

-3.81 

0.  18 

4.71 

3ii  76 

147 

-2.  92 

0.21 

4.  62 

3.97 

157 

Li 

11 


-3.  80 
-3.48 


0.20 
0.  23 


4.80 

4.72 


3.78 

3.89 


119 

131 


and  are  intended  to  be  uaed  only  ae  a vieual  aid  in  separating  earthquakes 
' from  presumed  nuclear  explosions.  Little  emphasis  should  be  placed  on 
their  particular  slopes  and  intercepts. 

Copsidering  the  individual  events  of  Table  VI- 1 for  which  sur- 
face wave  magnitudes  could  be  measured,  the  following  points  should  be 
noted  from  Figures  VI-2  to  VI- 14s 

• Event  839  The  surface  wave  magnitudes  for  this  event 

> 

fall  well  within  the  earthquake  population. 
Therefore,  based  on  this  discriminant,  this 
event  is  classified  as  an  earthquake. 

e Event  958  The  surface  wave  magnitudes  for  this  event 

fall  at  the  lower  edge  of  the  earthquake  pop- 
ulation. Therefore,  based  on  this  discriminant, 
this  event  cannot  be  classified  as  an  earth- 
quake or  explosion. 

• Event  1368  The  surface  wave  magnitudes  for  this  event 

show  excellent  separation  from  the  earthquake 
population.  Therefore,  based  on  this  discrim- 
inant, this  event  is  classified  as  an  explosion. 

e Event  1558  The  surface  wave  magnitudes  of  this  event 

show  excellent  separation  from  the  earthquake 
population.  Therefore,  based  on  this  discrim- 
inant, this  event  is  classified  as  an  explosion. 

e Event  1624  The  surface  wave  magnitudes  for  this  event 

fall  well  within  the  earthquake  population. 
Therefore,  based  on  this  discriminant,  this 
event  is  classified  as  an  earthquake. 


VI-21 


SECTION  vn 

COMPARISON  OF  SINGLE-SITE  STATION 
AND  ARRAY  EVALUATION  RESULTS 

The  goal  of  this  section  is  to  determine  the  differences  in 
performance  between  the  seven-element  Iranian  Long-Period  Array  and  a 
single-site  station.  The  single-site  station  chosen  for  this  comparison  is 
Mashhad  (MAIO)  of  the  Seismic  Research  Observatory  Network,  located 
approximately  700  km  from  ILPA.  This  station  was  selected  because  it 
is  closest  to  ILPA  of  the  available  single- site  stations  and  has  the  same 
sensor  system  as  the  individual  remote  sites  of  ILPA.  All  data  for  MAIO 
are  taken  from  the  current  report  on  the  evaluation  of  the  Seismic  Research 
Observatory  stations  (Strauss,  1977). 

If  the  problems  which  have  shut  down  ILPA  in  the  past  as  de- 
scribed in  Section  II  do  not  recur,  the  reliability  of  the  array  as  a whole  is 
assessed  as  0.  8 (where  a reliability  factor  of  1.  0 would  be  perfect. ) At 
MAIO,  this  reliability  factor  is  assessed  at  0.  9.  The  difference  in  these 
reliability  estimates  can  be  ascribed  to  the  greater  complexity  of  instru- 
mentation at  ILPA,  resulting  in  a higher  probability  of  hardware  failure  and 
consequent  down-time  at  the  array. 

The  comparison  of  mixed  event  probabilities  is  presented  in 
Table  VU-1.  Note  that  while  at  MAIO  the  probability  of  an  event  being 
mixed  remains  fairly  constant  for  all  ranges  of  bodywave  magnitude  consid- 
ered, It  tends  to  drop  for  increasing  bodywave  magnitude  at  ILPA.  This 
difference  is  due  to  the  process  of  beamforming  array  data.  If  the  event 
causing  the  mixing  is  well  off  the  azimuth  of  the  event  under  analysis  and 
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if  the  amplitude  a of  the  two  events  as  recorded  at  the  array  are  roughly  the 
same,  beamforming  will  suppress  the  off-azimuth  signal  and  reveal  the  sig- 
nal under  analysis.  This  change  from  mixed  event  on  single- site  to  detected 
event  on  beam  occurred  for  events  processed  at  the  higher  values  of  body- 
wave  magnitude. 

Table  VXX-2  presents  the  comparison  of  RMS  noise  levels 
measured  in  the  0.  023  - 0.  059  Hs  passband  at  MAIO  and  ILPA.  Note  first 
that  the  RMS  noise  at  MAIO  is  slightly  less  than  the  RMS  noise  measured  at 
the  ILPA  reference  site.  Beamforming  the  noise  recorded  at  ILPA  lowers 
the  ILPA  noise  levels  by  4.  7 dB  for  the  vertical  component,  3.4  dB  for  the 
north  component,  and  2.  0 dB  for  the  east  component  relative  to  the  MAIO 
RMS  noise  levels.  For  the  two  highest  components,  this  is  an  average  noise 
suppression  of  4.  0 dB  relative  to  MAIO.  (This  averaging  is  justified  by  the 
consideration  that  an  event  is  considered  to  be  detected  if  it  is  detected 
on  two  of  the  three  components  and  that  it  is  most  likely  to  be  detected  on 
the  components  with  highest  noise  suppression. ) From  this  on#  can  predict 
that  the  detection  capability  of  ILPA  relative  to  MAIO  (as  measured  by  the 


50  percent  detection  threshold)  should  be  about  0.  20  m^  units  lower. 

Table  VII- 3 presents  the  comparison  of  MAIO  and  ILPA  de- 
tection capability.  The  values  of  mean  epicentral  distance  ('MEAN  DELTA') 
are  presented  as  an  indication  that  the  MAIO  and  ILPA  data  bases  from  which 
the  detection  capability  estimates  were  determined  had  overall  essentially  the 
same  epicentral  distances.  Therefore,  the  comparison  is  not  obscured  by 
differences  in  detection  capability  due  to  differences  in  epicenter- station 
separations.  The  terms  'ALL  REGIONS  A'  and  'ALL  REGIONS  B<  are  as 
previously  defined  in  Section  V,  where  'ALL  REGIONS  A'  represents  the  ab- 
solute detection  capability  estimate  and  'ALL  REGIONS  B'  represents  the  con- 
ditional detection  capability  estimate. 
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TABLE  VH-2 


COMPARISON  OF  RMS  NOISE  LEVELS 
(ALL  VALUES  IN  mfi) 


Station 

Vertical 

North 

East 

Mean 

S.D.  * 

Mean 

S.D.  * 

Mean 

S.D.* 

MAIO 

8.20 

4.37 

7.89 

3.45 

8. 19 

4.  00 

ILPA  Reference  Site 

9.58 

2.19 

9.  34 

2.  56 

9.86 

2.  98 

ILPA  BEAM 

4.77 

. 

1.44 

5.36 

1.82 

6.  51 



2.31 

* S.D.  * Standard  Deviation 


Note:  RMS  noise  computed  in  0.  023  - 0.  059  Hz  paiBband 
(no  correction  for  instrument  response) 
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TABLE  VII- 3 

COMPARISON  OF  ILPA  AND  MAlO  DETECTION  CAPABILITY 


Type  of 
Detection 
Statistic* 

MAIO 

ILPA 

®b90 

Mean  Delta 
(degrees) 

"Vo 

Mean  Delta 
(degrees) 

All  Regions  A 

All  Regions  B 

4155 

381  5 

4.42 

36.5 

4. 13 

3.94 

VH-5 


The  data  of  Table  VII- 3 show  that  the  ILPA  detection  capability 
is  affected  slightly  more  than  the  MAIO  detection  capability  by  the  aggregate  of 
mixed  events,  malfunctions,  and  system  failures  causing  no  data  to  be  recorded. 
(Note  that  ILPA  is  evaluated  only  for  those  times  when  it  was  nominally  oper- 
ational. ) For  ILPA,  the  absolute  and  conditional  detection  capability  estimates 
differ  by  0.48  m^  units,  while  for  MAIO  the  difference  is  0,42  units.  Con- 
sidering the  mixed  event  probabilities  of  Table  VII- 1,  this  difference  must  be 
a direct  result  of  the  lower  reliability  of  ILPA. 

A value  of  0.  20  units  was  predicted  from  the  relative  ILPA 
and  MAIO  RMS  noise  levels  as  the  improvement  in  the  ILPA  detection  capabil- 
ity relative  to  the  MAIO  detection  capability.  Since  thiB  does  not  take  into 
account  the  effects  of  mixed  events,  malfunctions,  and  system  failures  result- 
ing in  no  data  being  recorded,  it  represents  an  improvement  in  the  conditional 
detection  capability  estimates.  From  Table  VII-3,  the  ILPA  conditional  de- 
tection capability  estimate  ('ALL  REGIONS  B1)  is  0.  19  m^  units  lower  than 
that  of  MAIO,  which  agrees  quite  well  with  the  predicted  improvement. 


SECTION  VIII 
CONCLUSIONS 


This  section  summarizes  the  results  of  the  continued  eval- 
uation of  the  Iranian  Long-Period  Array  (ILPA)  and  presents  plans  for  nec- 
essary future  work.  The  major  conclusions  are: 

A.  DATA  QUALITY 

• In  general,  the  data  quality  is  fairly  good.  Of  680  events 
examined,  6.  6%  were  lost  due  to  unreadable  data,  6.  2%  were 
lost  due  to  gaps  In  the  recorded  data,  and  4. 1%  were  lost  due 
to  uncorrectable  system  malfunctions. 

• Estimates  of  array  reliability  were  based  on  estimates  of 
array  up  time,  frequency  of  occurrence  of  recorded  and  read- 
able data,  and  frequency  of  occurrance  of  system  malfunctions. 
The  worst  case  estimate  of  ILPA  reliability  is  0.  54  and  the 
best  case  estimate  is  0.  79. 

e The  probability  of  an  event  being  mixed  was  estimated  at  approx- 
imately 0.  25  for  a range  of  m^  values  from  3.  1 t0  5.  0. 

B.  NOISE  ANALYSIS 

• RMS  noise  amplitudes  do  not  vary  greatly  from  site  to  site. 

• Noise  suppression  due  to  beamforming  is  greater  on  the  vertical 
component  than  on  either  of  the  horizontal  components. 

e For  all  sites  except  site  1,  the  RMS  noise  amplitudes  for  periods 
above  25  seconds  are  higher  on  the  horizontal  components  than 
on  the  vertical  component. 
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At  periods  above  25  seconds,  beamforming  lowers  the  vertical 
component  RMS  amplitude  spectrum  much  more  than  the  hori- 
zontal spectra. 

The  above  points  indicate  that  use  of  a 0.  033  - 0.  050  Hz 
passband  filter  in  place  of  the  0.  023  - 0,  059  Hz  passband 
filter  previously  used  may  significantly  enhance  the  signal-to- 
noise  ratio  of  the  ILPA  data. 

Multichannel  noise  coherencies  for  the  7- element  1LPA 
array  are  quite  similar  in  level  as  a function  of  frequency  to 
those  determined  from  the  Inner- ring  sites  of  the  Norwegian 
Seismic  Array. 

The  level  of  multichannel  coherence  in  the  0.  023-0.  059  Hz 
filter  passband  is  high  enough  that  multichannel  filtering  may 
be  effective. 

There  is  very  little  0,  01172  Hz  propagating  noise  arriving 
at  ILPA. 

The  majority  of  0.  04297  Hz  and  0.  05859  Hz  propagating 
noise  with  signal  phase  velocities  have  arrival  azimuths  dir- 
ected toward  the  south,  away  from  the  general  seismic  area 
of  interest. 


SIGNAL  ANALYSIS 


The  highest  gain  ip  signal-to-noise  ratio  due  to  beamforming 
was  6.  9 dB  on  the  vertical  component.  Gains  on  the  trans- 
verse and  radial  components  were  4.4  and  4.  0 dB,  respectively. 

The  LQ-transverse  component  yielded  larger  correlation 
coefficients  than  either  the  LR-vertical  or  LR-radial  components. 
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The  LR- radial  component  produced  the  lowest  correlation 
coefficients  on  average.  Sites  3 and  6 gave  the  highest  corre- 
lation coefficients,  while  site  4 yielded  the  lowest  correlation 
coefficients. 

DETECTION  CAPABILITY 

e The  absolute  50  percent  detection  capability  estimate  for  ILPA 

beam  data  using  NORSAR  m^  values  is  at  m^  = 4.42  for  Eurasian 
events.  The  absolute  detection  capability  estimate  was  computed 
by  including  all  mixed  events,  events  for  which  no  data  were  avail- 
able and  events  containing  malfunctions  as  non-detectione. 

e The  conditional  50  percent  detection  capability  estimate  for  ILPA 
beam  data  using  NORSAR  values  1b  at  m^  * 3.  94  for  Eurasian 
events.  The  conditional  detection  capability  estimate  was  com- 
puted by  excluding  all  mixed  events,  events  for  which  no  data 
were  available  and  events  containing  malfunctions  from  the  detec- 
tion statistics. 

e Mixed  events,  events  for  which  no  data  were  available  and  mal- 

functions raised  the  50  percent  detection  threshold  for  the  ref- 
erence site  and  beam  by  0.41  and  0.48  m^  unitB,  respectively. 

DISCRIMINATION 

• Events  839  and  1624  show  surface  wave  magnitudes  indicating 
that  they  are  earthquakes, 

• Event  958  could  not  be  classified  as  either  earthquake  or  ex- 
plosion. 

• Events  1368  and  1558  show  surface  wave  magnitudes  indicating 
that  they  are  nuclear  explosions. 
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F.  COMPARISON  WITH  A SINGLE-SITE  STATION 

• ILPA  reliability  is  estimated  at  0.  8,  in  comparison  with  Q.  9 

at  MAIO.  This  reflects  the  greater  system  complexity  of  ILPA. 

• Due  to  the  beamforming  process,  the  probability  of  an  event 
being  mixed  decreases  with  increasing  bodywave  magnitude  at 
ILPA  while  remaining  fairly  constant  at  MAIO. 

• While  the  RMS  noise  levels  at  ILPA  are  slightly  greater  than 
those  at  MAIO,  beamforming  produces  an  average  noise  suppres- 
sion of  4.  0 dB  at  ILPA  relative  to  MAIO. 

• The  conditional  detection  capability  estimate  at  ILPA  is  approx- 
imately 0,  2 m,  units  lower  than  that  of  MAIO. 

b 

G.  FUTURE  WORK 

The  following  items  should  be  investigated  to  complete 
the  evaluation  of  the  Iranian  Long- Period  Arrays 

• Investigate  the  long-period  noise  field  in  more  detail.  Points 
covered  should  include  RMS  noise  levels  in  different  passbands, 
reasons  for  increased  horizontal  component  noise  at  periods  be- 
yond 30  seconds,  and  frequency-wavenumber  spectra  for 

more  noise  samples  to  better  define  source  azimuths  of  prop- 
agating noise. 

• Refine  regionalized  detection  capability  estimates.  This  will 
require  processing  and  analyzing  approximately  500  more  events. 

• Process  more  presumed  nuclear  explosions  to  better  estimate 
the  ILPA  discrimination  capability. 
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• Process  a suite  of  detected  and  non-detected  events  with 
the  0.  033  - 0.  050  Hz  passband  to  determine  the  effect  of 
removing  the  elevated  noise  levels  beyond  30  seconds  on 
detection  capability. 

e Completely  evaluate  the  single  recorded  component  of  short- 

<i'  j | ' period  data  in  termB  of  noise  characteristics,  detection 

capability,  and  discrimination  capability. 
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APPENDIX  A 
THE  DATA  BASE 


This  appendix  presents  the  parameters  describing  each  event 
in  the  data  base.  The  column  headed  'EVNO*  gives  the  number  assigned  to 
each  event.  (These  event  numbers  are  the  same  as  those  used  in  the  Seismic 
Research  Observatory  evaluation. ) The  column  headed  'DATE'  gives  the 
month,  day,  and  year  of  occurrence  oi  the  event.  The  column  headed  'TIME' 
gives  the  origin  time  of  the  event.  The  columns  headed  'LAT.  1 and  'LONG.  1 
give  the  latitude  and  longitude  of  the  event  epicenter,  where  a positive  value 
Indicates  north  latitude  or  east  longitude  (as  appropriate)  and  a negative  value 
indicates  south  latitude  or  west  longitude.  The  column  headed  'MB'  gives 
the  body-wave  magnitude  of  the  event.  The  column  headed  'Q'  gives  the 
NORSAR  quality  rating  of  the  event  parameters,  where 

1 = good  to  excellent 

2 * fair  to  good 

3 a poor  to  fair. 

A zero  in  this  column  indicates  the  event  parameters  came  from  the  NEIS 
event  bulletin.  The  column  headed  'LOCATION'  giveB  the  general  area  in 
which  the  event  occurred.  Finally,  the  column  headed  'SUBREG'  gives  the 
sub-region  number  of  the  event  as  d<  fined  by  Flinn  and  Engdahl  (Flinn  and 
Engdahl,  1965). 
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